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ABSTRACT 


The current antenna used in the Marine Radio AN/MRC-142 is a focal point fed 
reflector with a first sidelobe of -10 dB. A new antenna design is required with a sidelobe 
reduction of another 10 dB, and a minimum impact on the system upgrade cost. An 
improved feed consisting of a fed dipole and disk in the vicinity of a parasitic ring was 
designed for use with the current reflector. The disk acts as a ground plane for the dipole, 
with the ring narrowing the beamwidth, resulting in a stronger edge taper on the reflector. 

An existing computer code based on the the method of moments was used to 
compute the gain and the radiation pattern. Circular parabolic reflectors of 36 inches and 
48 inches were analyzed. The results show that the 36 inch reflector cannot satisfy the 
conditions over the entire frequency band (1350 to 1850 MHz) simultaneously. The 48 
inch reflector came very close to satisfying all the conditions, although the sidelobe level 


is off by approximately 2 dB at wide angles in the H-plane. 


1 


TABLE OF CONTENTS 


A. THE E-FIELD INTEGRAL EQUATION AND THE METHOD OF 


MOMENT 8 occa cas eelek ee a ct 


l; Baseline Feed) Desion 2.225 1. eee 


2. Disk Size And Location... .. eee 


o0: me mw 


3. Ring Size And Position ........... - re ‘ a ee 


Pome NOI Oedel abasic IOIDOlen ess... s,s ss. ee On 
pmberocussing Nie Need 5. 2 oe oes ee vaste vvvesussestess oe 35 
Pee NOM MOINES SIGN Se ce ee bec eee webs utes eenme des a5 
mesicle Prequency (1550 Mhz) 2.2... 02. sek ee see eee ee. nee 3)5 

2. Dual Frequency (1350 And 1850 Mhz) ....................... ee 38 

IV. CONCLUSIONS AND RECOMMENDATIONS ..........0..0000000000002.0..: 45 
APPENDIX A. ANTENNA SPECIFICATIONS ...........000000.0.0 00000 0c. 47 
Piel B. PROGRAM CODE .............. 0.0.0.0 cece ccc cece ee ee eee 49 
TSE EIEN ES oc. ccc ce eee else eee yd cee eee e seen seed dt a neee ees 68 
Reet OISTRIBUTION LIST ...... 2.22. cece 69 


ACKNOWLEDGMENT 


Completion of this research effort would not have been possible without the help of 
Dr. David C. Jenn and Dr. Rama Janaswamy. Their guidance has extended beyond the 
scope of this work to my entire stay at the Naval Postgraduate School, and I am deeply 
appreciative. Additionally, I would like to acknowledge the consideration of Carmen 


Cibeira, whose understanding and support was a source of strength. 


v1 


I. INTRODUCTION 


A. BACKGROUND 

An antenna 1s the transducer that intercepts propagating electromagnetic waves and 
directs them into guided waves in a transmission line. In the case of a transmitting 
antenna, guided electromagnetic energy in a transmission line is converted into radiated 
electromagnetic waves. In short, the antenna couples the system to the environment [Ref. 
1]. 

Almost all radar, communication, and electronic warfare (EW) systems require 
some form of antenna. Even for highly directive antennas, all of the energy is not 
concentrated into the main beam; there is always radiation in directions other than the 
intended one. These sidelobes, as they are called, are of particular importance in 
electronic warfare, since they represent a major limitation to the system's effectiveness. 
Jamming signals or false targets received through a sidelobe can cause confusion and 
failure of a mission more easily than false signals detected by the main beam. The 
reduction and control of the position and magnitude of sidelobes is extremely important 
when attempting to localize target azimuths 1n a hostile electronic environment. 

Antenna sidelobe control is one of the most basic methods of defeating jammers 
and avoiding detection. Low sidelobes also reduce multipath signals, which can cause 


signal fading in communications systems and tracking errors in radars. The methods of 


sidelobe suppression for passive antennas are well known, and the technology is low risk. 
However, low sidelobes are accompanied by a reduction in gain, and generally require a 
more complicated antenna design. Consequently, there are tradeoffs in the antenna design 


that impact the overall system performance and cost. 


B. OBJECTIVE 

This thesis investigates the feasibility of reducing antenna sidelobes for the Marine 
Radio MRC-142. This system provides medium-range ultra high frequency (UHF) 
line-of-sight multi-channel radio communications over a distance of up to 35 miles. The 
AN/MRC-142 radio terminal set is installed in the M-998 variant of the High Mobility 
Multipurpose Wheeled Vehicle (HMMWV) and when deployed appears as in Figure 1. 

This unit can be configured as a repeater to extend the communications range 
between two other systems or to overcome terrain conditions that prevent a direct 
line-of-sight path between the other two pieces of equipment. The AN/MRC-142 
interfaces with a distant end AN/MRC-142 over a 1350 to 1850 MHz UHF radio link 
using parabolic antennas that are mounted on masts which can be extended up to 50 feet. 
The current design uses a parabolic reflector with a focal-point feed and has a first 
sidelobe requirement of -10 dB relative to the mainbeam. A reduction of another 10 dB is 


desired for the new design, with a minimum impact on the system upgrade cost. 





50 Foot Telescoping Antenna 





Rack Mounted Equipment (RME) 





M-998 HMMWV 


Figure 1. Deployed AN/MRC-142 Radio Terminal Set [Ref. 2]. 


C. METHOD OF SOLUTION 
Possible options for the low sidelobe antenna include: 
]. an array antenna; 
2. acompletely redesigned reflector and feed, or; 


3. an improved feed used 1n conjunction with the current reflector. 


The complete antenna specifications are given in the Appendix. Options | and 2 could 


satisfy the specifications in their entirety, but the cost of design and development would 


be significant. The third option is the one pursued here. Even though all of the antenna 


specifications are not achieved simultaneously with the new feed, there is significant 


improvement in all areas. Most importantly, the cost of the new feed is minimal because 


of its simplicity. 


The proposed feed design consists of a fed dipole in the vicinity of a parasitic ring 
and disk as shown in Figure 2 [Ref. 3]. Essentially the disk acts as a ground plane for the 
dipole and the ring narrows the beamwidth, resulting in a stronger edge taper on the 
reflector. The dimensions are chosen to minimize the sidelobes. This "tuning" approach 
is narrow band and, as discussed in Chapter III, does not allow an optimum feed design 


for both the transmit and receive frequencies. 


DIPOLE 


DISK 


Figure 2. Dipole-Disk Antenna with Ring Used as a Feed in Paraboidal Reflector. 


D. SCOPE OF THESIS 
In Chapter II the method of analysis is discussed in detail. All of the gain and 


pattern data presented are generated using the method of moments. Chapter III discusses 


the tradeoffs involved in the antenna design. As mentioned earlier, the greatest challenge 
was selecting the geometrical parameters to provide acceptable performance at both 
frequencies (1350 MHz and 1850 MHz). Finally, Chapter IV presents a summary of the 


work performed along with some conclusions. 


Il. THEORETICAL BACKGROUND 


A. E-FIELD INTEGRAL EQUATION AND THE METHOD OF MOMENTS 

To determine the pattern and gain of a reflector antenna, it is first necessary to find 
the current on the reflector surfaces. This 1s accomplished by solving an integral or 
differential equation for the current which is derived from Maxwell's equations and the 
boundary conditions. For wire and surface geometries of arbitrary complexity, the 
integral or differential equation must be solved numerically on a high speed digital 
computer. The method of moments (MM) ts such a technique which reduces the E-field 
integral equation (EFIE) to a matrix equation which 1s then solved on a computer [Ref. 4]. 
The following paragraphs summarize the method of moments technique. Phaser notation 
is used throughout and e’ time convention is assumed and suppressed. 

If the current on a body is known, then the scattered (or radiated) field can be 


determined from [Ref. 5] 
E=-Vo-joA =-joA-janV(V eA), (1) 


where ® = 27f is the radian frequency, and and € are the permeability and permittivity 


of free space. A is the magnetic vector potential, which can be represented as 


— t = / / 7 e~JkR f 
A==[,Js(« .y ,2 )-ds (2) 


where J. is the surface current and S refers to the surface of the conductor. 


The electric field integral equation is based on the boundary condition that the total 
tangential electric field on a perfectly electric conducting (PEC) surface of an antenna or 


scatterer is zero. This can be expressed as 


Ei¢ =T,) = E\@ =1;) +ES@ =Ts)=0 onS (3) 


where S is the conducting surface of the antenna or scatterer, as stated before, andr =r, 


is the distance from the origin to any point on the surface of the antenna or scatterer. The 
subscript "t" indicates tangential components. 


Referring to Figure 3, Equations (1) and (2) can be combined and expressed as 





Figure 3. Coordinate Systems for Computing Radiation Fields. 


Es(r) = -j71B7J J, Js@ GET )ds’ + VI J, V! © J@)GE,T ds’), (4) 


where 


eJkR 


Cer) == (5) 
and 

R= |r -T’]. (6) 
Now, applying the boundary condition of Equation (3) gives 


E°(r) =E%(T,;)= [Bf J, JF G(7,7ds/ + VJ J, V’ eJ.(r’)G(T,1)ds’ |) ae 
This is an integral equation with J, as the unknown. Equation (7) is referred to as the 


electric-field integral equation. In Equation (7), V and V’ are the gradients with respect to 


the observation (unprimed) and source (primed) coordinates, respectively, and 
/ 
Gr rs) 
is the Green's function for a three dimensional radiator or scatterer. 


For convenience, an operator L can be defined which will operate on J 
LU.) =jousJ.J.G; ds’ +5V| Vif, J.G; ds’ +J JJ, VG; ds’ (8) 


where unknown current J, depends only on the primed quantities. G is a scalar quantity 


that depends on both the primed and unprimed quantities. It 1s desirable to have all the 
differential operations in terms of the primed variables of integration. Equation (8) can 


therefore be written as 


LG.) =jonsl,IsG; ds’ +EV| ff, G; VeTsds’+JJ,J,°V; G; ds’, (9) 


where the vector identity 
V(J, ¢ VG) = (J, ¢ V)VG+ (VGe V)J, +3, x(V x VG)+VGx(V xJ,). (10) 


has been used. 


The gradient of G in terms of the two sets of coordinates is related by 


—)kT a 





VG =-(jk+7)5 


Ate 


=-V’/G. (11) 


Using the surface divergence theorem developed by Mautz and Harrington [Ref. 5], it can 


be shown that 

[{. V/ «© (I,G)ds’ = 0. (12) 
Finally, Equation (7) can be written as 

LO.) =J {Lion .G-EVE(V! © J.)G] Ids’. (13) 
Now the electric field integral equation can be represented as 

Ba.) =Ld,). (14) 


Equation (14) can be solved numerically using the method of moments. The current 
on the antennas is expanded in terms of a finite series of known basis functions with 
unknown amplitudes. The antenna under consideration consists of surfaces and thin 
wires, so basis functions must be defined on each. On a surface two orthogonal 
components are required. The coordinate system for the reflector 1s shown in Figure 4 


and the following basis functions defined on the surface [Ref. 6]: 


SEGMENTED REFLECTOR GENERATING CURVE 


Pitt 





| 2 
2 int Zi AXIS OF SYMMETRY 


Figure 4. A Reflector Geometry for Method of Moments. 


—_— aT, _: 
Jt, =t-ein , ie O Sete. o «eae 
Hz iOe) 


Jo = Geen, m=O 41... 


(15) 


T(t) 1s the triangle function (which extends over two segments, 1 andi + 1) and P(t) 


is the pulse function. A point on the surface of the antenna 1s specified by the coordinates 


(t, >), where t is an arclength variable along the reflector generating curve. The distance 


of a point from the axis of symmetry (z-axis) 1s given byp and N, is the number of 


surface generating points. A discussion leading to the choice of these basis functions 1s 


given by Mautz and Harrington [Ref. 5]. 


For thin wires the radius "a" is much less than a wavelength and the current will 


only have an axial component. Thus the operator V becomes t(O/Ct), where t 1s the arc 


length along the wire. Furthermore, the current is constant around the wire and SJ,ds’ 
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reduces to 27aJ,dt’, as shown in Figure 5. Thus appropriate basis functions on the wire 


are [Ref. 7]: 


a See ee 


Figure 5. Thin Wire Approximation. 


ii i=1,2....N,. (16) 


oma ? 





In this case t is a unit vector along the wire axis and N,. is the number of wire generating 


points. The total current is the sum over all expansion functions 


= Nw-2 0 oo Ns-2 a Ns-1 6 =4 
ie= > lLJnt > 2 eae Ds JgnJ til (17) 
m=! N=—-C} p=1 q=l 


The coefficients I are the unknowns which need to be determined. 
From the linearity of L, Equation (14) can be combined with Equation (17) as 


follows 


an N = 
Ee » LLG). (18) 


where N, is the total number of basis functions of both types. N, equations need to be 
provided to enable a solution for the unknown constants I,. Galerkin's method is used to 


reduce Equation (18) to a matrix equation, and weighting functions are defined such that 


Wi = Ts (Ref. 8, 9]. Each side is multiplied by each of the basis functions and then each 


side integrated. The result 1s N, equations of the form 


[J Ee Weds = [JJ Ws eli Li) +J J, Wee 12LG2) +... +I J, Ws ¢ InLGn)] ds 
ieee (19) 


which can be written as 
= = Ss Np = 
JJ, Bie Wids=J J, [Wee = 1,L( ;)ds]. (20) 
Interchanging summation and integration, 
260.2 Ny = 2 
ie E'e Wids = & Tif] Wx e Lids. (21) 
Equation (21) can be cast in the form 
Nb 
Wie Ze (22) 
i=] 
for k = 1, 2, ...,N,, where the impedance elements are given by 
Zu =Sf,dsJf,|jouWs Ti-Z(V © FV © Wy) [Gas | (23) 
The excitation elements are defined by 


Vi =JJ. Wee E'ds, (24) 
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In matrix notation Equation (22) becomes 
LV} = (2), (25) 
and the unknown coefficients are determined by solving the matrix equation 


 =[Z}"[V]. (26) 


Submatrices in Equation (25) associated with the t and b surface components can be 


identified as follows 


P=) (ZoZal || 1 
| [Va] i: | ZZ | | ie | em 


The impedance matrix will have the following block structure when all of the 


surface and wire basis functions are considered 


SS SW S S 
ge ee ee lew, Wien, 
SS sw S s 
SW S s 
Zy1 ) Poet ty Ya 
“titi WSo 8 Ws ww oi | W W | 
bee Ze ZZ aL A Ve | 


The superscripts s and w refer to surface and wire, respectively. Although n in 
Equation (17) can take on values between + %, it has been truncated at+ M in the above 
matrix. In practice, when the feed 1s on the reflector axis of symmetry, a converged 


solution can be obtained with M = 1. If the pattern is calculated for the antenna 
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transmitting, the excitation vector will have only one non-zero value corresponding to the 
dipole feed point. 

The unknown expansion coefficients are determined by solving the matrix in 
Equation (27) for [I]. To compute the electric field, the coefficients are used in Equation 


(17), which in turn is used in Equations (1) and (2). Define a measurement vector, 


— — 


Rm =Jf Ere Jsds, (29) 


where E, denotes a unit radiated plane wave weighted by the current on S. In general Ui, 


will give rise to both 8 and@ components in the far field. The radiated 8 component is of 


the form 
E, es Beniker — He 7iklsutyvtzw] (30) 


where u = sin@cosd, v = sin@sinsd, and w = cos0. By substituting Equation (30) into 


Equation (29), there results 


R° a2 ae e“iklsutyvtzw] (6 eJ.)ds . é 1) 


The total 8 component of the electric field is the superposition of all surface current 


contributions 
-jkn on 
Es» = a a Ro tee (32) 
The analysis holds true for the @ component. The measurement elements are 


Rh = ff ether“ (ho J)ds , (33) 


and the electric field is then defined by 


Ea ery Ree (34) 


4nr 


B. GAIN CALCULATION 
The degree to which any particular antenna pattern is concentrated into a beam is 


known as directivity which 1s defined by 


D=—= (35) 


Pave ; 


where P,,,, 18 the maximum power flux radiated, and P,,. is the average power flux 
radiated. Since the average power flux radiated is equal to the total power radiated (P.,,,) 


divided by 47 steradians, 


Pinax 
2 i (3 6) 


Directivity, as defined above, is a function of the antenna pattern only, and does not take 
into account the antenna efficiency or losses. In order to compare the relative 
performance of antennas in any electronic-warfare system, it 1s necessary to consider 


antenna efficiency. The pertinent quantity 1s gain G which 1s defined as 


G=KD, (37) 


le 


where K is the efficiency factor (O<SK< 1). Thus it is seen that for a 100% efficient 
antenna (K = 1), the gain would be equal to the directivity. For all practical antennas K < 


1. Gain is frequently expressed as a decibel ratio 
G, db = 10log,,G. (38) 


A very important relationship exists between gain, physical antenna dimensions, 


and electrical wavelength, expressed by 
G= 4n®>, (39) 


where A is the "effective" area of the antenna andA 1s the wavelength. The effective area 
of the antenna and the wavelength must be expressed in the same units (meters, inches, 
feet, etc.). The total radiated power required in Equation (36) can be obtained by 


integrating the electric field over a sphere 


Pro = | [{|Eol? + |Es|7 }r2sin @dOd¢ (40) 
0 


0 


where E, and E, are known from Equations (32) and (34). 


Ili. LOW SIDELOBE ANTENNA DESIGN STUDY 


A. SUMMARY OF TRADEOFFS IN THE ANTENNA DESIGN 

The purpose of this thesis is to achieve a reduction in the antenna sidelobes for the 
Marine radio AN/MRC-142. The approach will be to upgrade the existing reflector in 
Figure 6 by implementing a feed design similar to that of Kildal and Skyttemyr [Ref. 3] 


shown in Figure 2 of Chapter I. 





Feed Assembly 





Antenna Assembly 


Offset Adapter 


Figure 6. AN/MRC-142 Antenna Assembly and Offset Adapter [Ref. 2.]. 


The size of the reflector aperture is 36 x 48 inches as shown in Figure 7. The 
projected aperture is not circular, but has two circular segments removed. The method of 
moments code that was used to compute the gain and radiation pattern can only handle 


rotationally symmetric reflectors. In the design of the antenna, circular parabolic 
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reflectors of 36 and 48 inches were analyzed. Additionally, the structure supports only 
vertical currents, whereas the method of moments code assumes radial and azimuthal 
currents (solid surface). The actual AN/MRC-142 performance will lie somewhere 
between the two. The pattern will resemble the 48 inch reflector in the plane of the wide 
dimension, and the 36 inch reflector in the plane of the narrow dimension. The effect on 
gain can be accurately estimated based on the area reduction from the 48 inch diameter 


circular antenna due to removal of the circular segments. 





Figure 7. Antenna Radius Determination. 


The reflector surface is a section of a parabaloid and 1s fed from the front. Four 


factors affect the efficiency of this type of antenna: 


1. Spillover efficiency; 


2. Illumination efficiency; 


3. Blockage by the feed; 
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4. Cable and line losses from the antenna input to the feed. 


Taking into consideration all these factors, most of the burden in maximizing the 
efficiency is placed upon the antenna feed. The dipole-disk antenna, which essentially 
consists of a halfwave dipole over a circular plane reflector, is widely used as feed for 
paraboloidal reflector antennas at VHF and UHF. A method of improving this feed 
consists in incorporating a circular conducting beam-forming ring (BFR) over the dipole 
in a plane parallel to the disk. The beam-forming ring compresses the H-plane pattern of 
the dipole-disk feed, whereas it has no significant effect on the E-plane pattern. When 
this antenna is used as feed for paraboloidal reflectors, it exhibits near-identical principal 
plane aperture illumination and reduces H-plane spillover lobes. 

The diameter of the disk should be as small as possible in order to reduce blockage 
loss. However a small disk and ring result 1s a broad feed radiation pattern which does 
not provide sufficient edge taper for a low sidelobe pattern. To obtain the best 
performance a parametric study was performed for the geometrical quantities shown in 
Figure 8 and listed in Table 1. The parameters and their effect in antenna radiation 
pattern and gain discussed in the following sections. All dimensions used by the 
computer code are in wavelengths. 

1. Reflector Size 

The size of the reflector is 36 x 48 inches as shown in Figure 7. Two different 
designs were obtained for each frequency: one for 36 inches and the other one for 48 


inches. At 1350 MHz, the wavelength is 8.75 inches and at 1850 MHz the wavelegth 1s 
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6.38 inches. Thus 36 inches corresponds to 4.1A and 5.6A, while 48 inches corresponds 
to 5.54 and 7.5 at these two frequencies. 

The focal-length-to-diameter ratio (f/D) of the existing reflector antenna could 
not be determined from AN/MRC-142 technical manuals. Typical values are in the range 
of 0.3 to 0.55. Smaller values are preferable to reduce spillover loss and increase edge 
taper. A value of 0.4 was chosen for f/D. This parameter could be varied slightly without 


significant modification of the antenna's gain and pattern. 







AW - WIRE RADIUS 


POCAL LENGTH 
FOD* SIAanETER OP PARABOLOID 


PAIN REFLECTOR 


Figure 8. Program Parameters 
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TABLE 1. DESIGN PARAMETERS 
DM = diameter of paraboloid 

DS = diameter of disk 

ZC = distance of dipole to ground plane 







FOD = focal length to diameter ratio of the parabolid 






ZR = distance of ring from fed dipole along the z axis 






RC = radius of ring 
SW = ring width (modeled as a thin strip) 
DLN = fed dipole length 
PLN = parasitic director dipole length 
DDIP = distance between DLN and PLN along the z axis 
AW = radius of dipoles 
HLGTH = parasitic (H plane) dipole length (not used) 











H plane dipole distance from feed along the y axis (not used) 





2. Gain Calculation 
Gain estimates for the reflector and feed were obtained by first integrating the 
total electric field to find the total radiated power P,,,. The result 1s then used in Equation 
[36] to determine the directivities shown in Table 2. These numbers include spillover and 
aperture efficiency because the total field is integrated (that is, feed plus scattered fields). 
Losses that are not accounted for and thus determine the value of K in Equation 
[37] include: 
1. Transmission line losses for the cable that feeds the dipole. Assuming 0.5 dB/ft 


and 3 ft of cable yields 1.5 dB of loss. 


ae 


2. Reflection loss from the reflector surface. This is due to the fact that the surface is 


a grid rather than a solid. A conservative estimate is 0.5 dB. 


3. Any impedence mismatch loss. A VSWR of 2:1 results is approximately 0.5 dB 


of loss. 


Based on the above numbers the total loss is about 2.5 dB. The corresponding 


gains are shown in Table 3. 


TABLE 2. CALCULATED DIRECTIVITIES. 


48 inches 36 inches 
1350 MHz 24.75 dB 22.19 dB 
1850 MHz 27 AAs 24.90 dB 


TABLE 3. CALCULATED GAINS WITH LOSS. 


48 inches 36 inches 
1350 MHz 22.25 dB 19.19 dB 
1850 MHz 23.94 dB 21.90 dB 





B. MODIFICATIONS TO ANTENNA DESIGN 
Each of the parameters of the feed and reflector were varied in small increments and 
their effect on the antenna pattern and gain was evaluated. The results are summarized 


below. 
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1. Baseline Feed Design 

The dimensions of the baseline feed design are listed in Table 4. The disk serves 
as a ground plane for the dipole. The phase center of the feed will vary with frequency 
and angle. If the disk is an effective ground plane the phase center of the dipole and its 
image will be approximately at the midpoint of the two. However in the following 
discussion the dipole is considered to be the center of the feed. 

Low sidelobes require a strong amplitude taper across the reflector, which in turn 
requires a narrow feed beamwidth. Generally a narrow feed beamwidth is associated with 
a feed aperture which results in increased blockage. Thus there is a tradeoff between 
sidelobe taper and blockage loss. 

The geometrical parameters in Table 4 were the starting point for the parametric 
study, which is summarized in the following sections. The baseline antenna performance 
is shown in Figure 10. The solid line is the E-plane (g=0°) pattern and the dashed line the 


H-plane (d=90°) pattern. 
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Figure 10. Baseline Antenna Pattern. 


TABLE 4. BASELINE FEED GEOMETRY 


Z distance or 
Variable Value coordinate’ 


4.1 LG = Zo 
0.7 Zi O88 
fed dipole length OS 27d = 0.25 


parasitic dipole length not part of baseline design 
fed dipole height above disk 05 Zc-Zd 
ring heigth above disk 0:25 Zc-Zr 


wire radius OO1A - 
i 





The fed dipole is located at Z = 0; positive Z is away from the reflector 
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2. Disk Size and Location 

The disk must be large enough to serve as an effective ground plane for the 
dipole. Generally the dipole height must be significantly less than the diameter of the 
disk. For the dipole and its image to add constructively in the direction of the reflector, 
the spacing should be approximately .254. The patterns for disk locations of Z = +.1A, 
and -.1A are shown in Figures 11 and 12. All other feed parameters are held constant 
(values given in Table 4). The reflector diameter is 4.14 and f/D = 0.5 in all cases. The 
high levels in the H plane patterns for 80° < 8 < 120° are due to feed spillover. A 
comparison of the figures shows that there 1s a tradeoff between spillover levels in the 
rear hemisphere (90 < 8 < 180) and sidelobe levels in the forward hemisphere (0 < 90°). 


The impact of changing the disk radius by .1A is illustrated in Figures 13 and 14. 





Relative Power, dB 





9, degrees 


Figure 11. Radiation Pattern with Disk Shift of +0.1A. 
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Figure 12. Radiation Pattern with Disk Shift of -0.1A. 
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Figure 13. Radiation Pattern with Disk Size .12. 
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Figure 14. Radiation Pattern with Disk Reduced -.1A. 
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3. Ring Size and Position 
The beam-forming ring (BFR) compresses the H-plane beamwidth but leaves the 
E-plane beamwidth essentially unaffected. Thus the ring size and location can be used to 
equalize the feed beamwidth. Patterns for two ring radii are shown in Figures 15 and 16. 
The ring acts somewhat like a director in a Yagi array when it 1s located in front 
of the dipole (towards the reflector) [Ref. 6]. When the ring is located behind the dipole it 
has little effect on the reflector scattered field. As illustrated in Figure 17 and 18, a shift 


of .1X results in about a | dB difference in the first sidelobe level. 
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Figure 15. Radiation Pattern with Ring Size of 0.6). 
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Figure 16. Radiation Pattern with Ring Size of 0.52. 
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Figure 17. Radiation Pattern with Ring at Z = -0.2A. 
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Figure 18. Radiation Pattern with Ring at Z = -0.12. 
4. Adding a Parasitic Dipole 
A parasitic dipole was added as a director in order to improve the directivity of 
the feed and therefore increase the reflector edge taper. This method is effective in 
narrowing the beamwidth of cavity-backed dipole antennas. In this case it was found that 
changing the size and position of the parasitic element has no significant effect on the 


performance of the antenna and does not offer any advantage. This is probably because 
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the parasitic dipole is not located within a cavity. Adding a second ring around the 
parasitic dipole could improve the performance. In an effort to keep the feed as simple as 
possible, the parasitic element was not considered further. 
5. Defocusing the Feed 

The phase center of the feed varies with frequency and angle, but is generally 
located between the fed dipole and disk. The antenna bandwidth can be increased by 
slightly displacing the feed from the reflector focus. Defocusing can cause a degradation 
at one frequency but an improvement at the other. Relative to an optimum pattern, 
defocusing will result in a spoiled beam; that 1s, null filling and change in beam slope. An 
example is shown in Figure 19. Defocusing also affects pattern levels in the rear 


hemisphere. 


C. OPTIMUM DESIGNS 
1. Single Frequency (1850 MHz) 
After making all the comparisons mentioned earlier, an optimum configuration 
was identified for the 36 inch parabolic reflector operating at 1850 MHz. The patterns are 


shown in Figure 20, and the configuration parameters shown 1n Table 5. 
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Figure 19. Radiation Pattern with Reflector Shifted 0.2A in the Positive Direction To 
Defocus the Feed. 
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Figure 20. Radiation Pattern of Optimum Single Frequency Design for 36 inch Reflector 
(The Dimensions are given in Table 6). 
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TABLE 5. OPTIMUM GEOMETRY FOR 36 INCH 
REFLECTOR AT LOW FREQUENCY 


BM ee 


This radiation pattern meets almost all the requirements listed in the Appendix 
A. The calculated total gain is G = 20.7775, which 1s sufficient, but the far out sidelobes 
are about 2 dB too high. Unfortunately this is a narrow band design and has poor 
performance at the lower frequency (1350 MHz). This design illustrates that achieving 
the specifications at a single frequency is relatively easy. However, broadbanding the 
antenna proves to be much more difficult. 

2. Dual Frequency (1350 and 1850 MHz) 

Starting from the configuration described in Table 4, attempts have been made to 
obtain a broadband design for the frequency range from 1350 to 1850 MHz. Figures 2] 
and 22 present the radiation patterns obtained for the 36 inch antenna which achieves 
most of the specifications given in the Appendix A. The gains are marginal at both 
frequencies and the beamwidth at the low frequency in too wide. Table 6 describes the 


configuration of this reflector. 


38 


a = ——. <== a, —— 





. ° . 
‘ ‘ ® e ° ° . . . 
eee rere rape Shee ese se gore ss eee ete wees e ee wee s oe cass ce ee creer ce cee eee ere eere eee teasers eree ere tse eseeeeves 
t \ s : ° e . « 
. . * ‘ ° . 
t . ° e . . ° 
. ® 


(PLANE 


Relative Power, dB 





0 20 40 60 80 100 120 140 160 180 


0, degrees 


Figure 21. Radiation Pattern of 36 inch Reflector for Lower Frequency (1350 MHz). 
Dimensions are given in Table 6. 
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Figure 22. Radiation Pattern of 36 inch Reflector for Higher Frequency (1850 MHz). 
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TABLE 6. OPTIMUM GEOMETRY FOR 36 INCH 
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For the 48 inch antenna, Figures 23 and 24 show the radiation patterns for the 


optimum configuration whose dimensions are listed in Table 7. 
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Figure 23. Radiation Pattern of 48 inch Reflector for Lower Frequency (1350 MHz). 


42 









 ceeseeiHePLANE! 48" REFLECTOR 
TEED pean oo tage ag ooo 





Me feos mse = cris ries tine eeaien ew sig #0 /6'sisicitisasgiase aiesiasalegitiss sie» 0246 ¢lesesie eestiaeeasd see ere see ees 
Ee -15 
ae 
of NN ete Ves aa ieee Wat eat Nis raise kes atWemn Al Sanaa sd ides Caslea erat ius 5 <dok ¢ 
zs 
2 Bb cag re ne 
ov 
12 
= -30 
O 
wa 





0 20 46 60 80 100 120 140 160 8 180 


0, degrees 


Figure 24. Radiation Pattern of 48 inch Reflector for Higher Frequency (1850 MHz). 
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It is worth noting that there is no feed configuration that can be used for both frequencies 





and simultaneously meet the specifications. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


The objective of this thesis was to find the optimum combination of the feed and 
reflector geometry over the frequency band 1350 to 1850 MHz. The current antenna size 
is 36 x 48 inches. Since the computer code only accepts circular reflector apertures, two 
Separate circular antennas were analyzed. One was 36 inches and the second was 48 
inches. The proposed feed design, which consists of a fed dipole in the vicinity of a 
parasitic ring and disk, was a good approach to achieve a significant improvement 
compared with the original geometry of the Marine radio AN/MRC-142. Because of its 
simplicity, the cost of the new feed represents a minimum impact on the system upgrade. 

Calculations for both sizes were made separately and, as it 1s shown in Chapter III, 
the 36 inch reflector could not satisfy the requirements at both frequencies 
simultaneously, but it gave the best result for the high frequency (1850 MHz). The results 
obtained with the 48 inch reflector come very close to satisfying all the requirements, 
although the sidelobe level (SLL) is off by approximately 2 dB at wide angles in the 
H-plane. 

An assumption with regard to the transmission line losses was made, which affect 
the total gain of the antenna. In the field calculations, solid surface reflection 
characteristics were assumed. In fact, the surface is a grid, and consequently the 


reflection coefficients for the 8 and @ components of the electric field will vary with the 
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angle of incidence. This effect will be small if the wire spacing is small compared to A. 
Achieving the required gain was not a problem for the 48 inch design, and small 
adjustments in the antenna feed geometry did not have a significant negative impact on 
the antenna performance. It was also determined that if there was an improvement due to 
an additional parasitic dipole, the increase in complexity did not merit its use. 

To summarize, a low sidelobe reflector antenna has been designed which fulfilled 
most of the requirements for the Marine radio. The dipole disk feed with beam-forming 
ring 1s electrically and mechanically simple. It can be supported by means of a stiff 
coaxial cable aligned along the axis of the paraboloid, and held by a center support which 
offers minimum blockage when compared with the usual strut supports, as was stated by 


Kildal and Skyttemyer [Ref. 3]. 
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APPENDIX A. ANTENNA SPECIFICATIONS 


Over the frequency range of 1350-1850 MHz, the antenna must have gains as 


described in Figure A. 


Gain 


-5 dB 


-195 dB 





Degrees 


10 20 30 


Gain (dBi 
AOC C“‘(C*(AX = Se note A 

At+10 TAR least 5 dB lower than A 
At+20, JA least 19.5 dB lower than A 
At+30to180 ‘IA least 27 dB lower thanA 


Shall satisfy the relationship 17.5 < A < 21 dB 


> 


Figure A. Gain Requirements. 
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TABLE A. ANTENNA/MAST REFERENCE DATA 
Antennna Lightweight, open grid, reflector 
1350 to 1850 MHz 


> 


Frequency Range 
Gain 


1350 MHz 
1850 MHz 
Breamwidth 
Horizontal 
7¥2 + 1° from boresight 
15 + 1° from boresight 
Sidelobes 
ntennna 
Vertical 
72 + 1°, -2° from boresight 
15 + 1°, -2° from boresight 
Sidelobes 
Polarization 
Weigth 17.5 Ibs 
Maximum Input Power 15W 
Maximum Output Power 2320 W Effective Radiated Power (ERP) 
45 in X 35 in. Depth: 16 in 


Mast Telescoping mast, consists of eleven 
carbon-fiber tubular sections 


18 dBi minimum 
20 dBi minimum 
+7.5° 


3 dB half power 
15 dB down 
10 dB min 


3 dB half power 

15 dB down 

10 dB min 

Horizontal or vertical polarization 


Reflector Dimensions 


Height Continously adjustable from 5 ft 8 in to 50 ft 


Continously adjustable from 5 ft 8 into 50 ft__ 
Length (for storage and transit) 
Maximum Wind (50 ft. erection) 60 mph (without ice) 

40 mph (with ¥ in ice loading) 


G) 


5 lbs 
15 minutes 


Weight 
2-Man Deployment Time 


Mast Accessories Time Stowage of mast accessories 


Total accessory bag weight: 44 lbs 


ae 
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APPENDIX B. 


Gain2t 


C PROGRAM mreint.f 
Cc 
C FAR FIELD PATTERN INTEGRATION FOR GAIN 
CALCULATIONS. 
C *** general wire/bor geometries can be handled °** 
Cc 
C >>>> INPUT FILE IS pcurrent - SPECIALIZED FOR 
PROGRAM mrc.f 
C >>>>HAS A GAIN NORMALIZATION LOOP 
INCORPORATED 
C >>>> WIRE DOES NOT HAVE TO LIE IN THE Z=0 
PLANE 
Cc 
DIMENSION 
RH(1000),ZH( 1000), AG(100),XG(100),XT(10),AT(10) 
DIMENSION 
Q(30),S(30) nw 1(2),nw2(2),pha(1000),xh(1000) yh(1000) 
COMPLEX 
RS(1800), ETW,EPW,U,CNT(1500),ET1,EP1,ET2,EP2 
COMPLEX 
CONJG,CMPLX,EXP1,EXP2,RW(2000),ET,EP.UC 
DATA pi/3.14159/ 
data nt/2/,xt(1),at(1)/.5773503,1./ 
c data 
nt/4/,xt(1),xt(10),at(1),at(10)/.33998104,.8611363115, 
c * .6521451548,.3478548451/ 
RAD=P1/180. 
Braco! 
W=(0),1.) 
UC=-U/4./P| 
OPEN(1,FILE='outgaus’) 
READ(1,°) NPHI 
NT2=NT/2 
DO 1 K=1,NT2 
K1=NT-K+1 
AT(K1)=AT(K) 
XT(K1)=XT(K) 
1 XT(K)=-XT(K) 
BO 3K=1,NPHI 
3 READ(1,*) XG(K),AG(K) 
OPEN(3 file='pcurrent’,STATUS='old’) 
READ(3,°) 
DM,AW,ZC,RC,DS,NP1,NP2,NBLOCK,NWIRES,FFAC 
READ(3,°*) 
NP,MT1,MP1,N.MT2,NSURF ,NROW,MHI,MODES 
write(6,°) 'NP1,NP2,N,MODES,NROW,NWIRES,NP=' 


write(6,*) NP1,NP2,N, MODES,NROW,NWIRES,NP 


READ(3,°) (xh(i),i=1,np) 
READ(3,") (yh(i),=1,np) 
READ(3,°) (zh(i),i=1,np) 
READ(3,”°) (rh(i),=1,np) 
READ(3,°) (pha(i),i=1,np) 
READ(3,°) (nw1(i),i=1,nwires) 
READ(3,*) (nw2(i),i=1,nwires) 


PROGRAM CODE* 


READ(3,°) (cnt(i),i=1,nrow) 
CLOSE(3) 
NS=NBLOCK’N 
OPEN(8,FILE='outint’) 
WRITE(8,8000) NP1,NP2,NBLOCK 
8000 FORMAT(//,10X,"** WIRE/BOR ANTENNA GAIN 
CALCULATION °"’, 
*//,10X,,NUMBER OF BOR POINTS (NP 1)=',i4,/,10x, 
“NUMBER OF DIPOLE POINTS (NP2)=',14,/,10X, 
“NUMBER OF MODE BLOCKS (2°M+1)='/I2) 
NDIVT=2 
NDIVP=2 
WRITE(6,°) ‘ENTER NDIVT AND NDIVP' 
READ(5,*) NDIVT,NDIVP 
$1=0. 
$2=180. 
Q1=0. 
Q2=180. 
WRITE(6,*) THETA START AND STOP ANGLES 
(DEGREES)’ 
READ(5,°) S1,S2 
WRITE(6,*) THETA START AND STOP ANGLES 
(DEGREES)’ 
READ(5,*) Q1,Q2 
C 
C DEFINE STANDARD BOR INDICES 
C 
MT1=NP1-2 
MP1=NP1-1 
MT2=NP2-2 
C NWIRES=1 SINCE DUMMY SEGS WILL HAVE ZERO 
CURRENT 
NWIRES=1 
NW1(1)=NP 141 
NW2(1)=NP 
C 
C GAIN NORMALIZATION LOOP (USE THETA=PHI=0) 
C 
ET 1=(0.,0.) 
ET2=(0.,0.) 
ETW=(0.,0.) 
CALL 
PLANEW(NWIRES,NW1,NW2,NP1,NP2,xh,yh,zh,0.,0.,RW) 
DO 610 L=1,MT2 
610 ETW=ETW+RW(L)°CNT(L+NS) 
EXP 1=(1.,0.) 
EXP2=CONJG(EXP 1) 
DO 611 M=1,MHI 
NM=M-1 
CALL PLANES(NM,NP1,0,NT,RH,ZH,XT,AT,O.,RS) 
NTOP1=MODES-NM 
NTOP2=NBLOCK-: (NTOP1+1) 
NS2=NTOP 1°N 
NS1=NTOP2*N 
DO 612 L=1,MT1 
ET 1=ET 1+RS(L)“CNT(L+NS1)°EXP1 
IF(NM.EQ.0) GO TO 612 
ET 1=ET 1+RS(L)°CNT(L+NS2)*EXP2 


* Although this is FORTRAN code, the column conventions are not observed for 


this listing so that space can be conserved. 


612 CONTINUE 
DO 613 L=1,MP1 
ET2=ET2+RS(L+MT1)°CNT(L+NS1+MT1)°EXP1 
IF(NM.EQ.0) GO TO 613 
ET2=ET2-RS(L+MT1)°CNT(L+NS2+MT 1)*EXP2 
613 CONTINUE 
611 CONTINUE 
ET=UC*(ET1+ET2+ETW) 
EMAX=CABS(ET) 
WRITE(6,*) "EMAX=’,EMAX 
WRITE(8,8008) NT,NPHI,NDIVP,NDIVT,S1,S2,Q1,Q2 
8008 FORMAT(//,10X,’NT=',13,/, 10x, 'NPHI=’,14,/,10X, 
* ‘NUMBER OF INTEGRATION INTERVALS IN 
PHI=’,i3,/,10x, 
*'NUMBER OF INTEGRATION INTERVALS IN 
THETA=',i3,/,10x, 
~ THETA LIMITS: START SWOP=—F7.2.2%:-7 27, 10x, 
*"PHULIMITS: START, SVOP= .F7.2.24,F7.2) 
Cc 
c theta integration intervals 
c 
$1=S1*RAD 
S$2=S2°RAD 
Q1=Q1°RAD 
Q2=Q2°RAD 
DS=(S2-S1)/FLOAT(NDIVT) 
DO 502 l=1,NDIVT+1 
S(l)=(I-1)°DS 
WRITE(6,”) *I,S()=",1,S(1) 
502 CONTINUE 
cS 
c phi integration intervals 
Cc 
DQ=(Q2-Q1)/FLOAT(NDIVP) 
DO 503 l=1,NDIVP+1 
Q(I=(I-1)°DQ 
WRITE(6,°*) *I,Q(l)=",1,Q(1) 
503 CONTINUE 
SUM=0. 
do 2000 JJ=1,NDIVP 
Filer. 
P2=(Q(JJ+1)+Q(JJ))/2. 
DO 2000 J=1,NPHI 
PHR=p1*xg(j)+p2 
DO 1000 Il=1,NDIVT 
T1=DS/2. 
T2=(S(Il+1)+S(II))/2. 
DO 1000 l=1,NPHI 
THR=T1°XG(l)+T2 
IF(THR.LT.Pl) GO TO 510 
THR=(BK-THR) 
PHR=PHR+PI 
510 CONTINUE 
STRA=SIN(THR)*AG(I)*ag(j) 
EFl1=(0-0,) 
EP1=(0.,0.) 
ET2=(0.,0.) 
EP2=(0.,0.) 
Elvv2(0.0)) 
EPW=(0.,0.) 
C 
C WIRE CONTRIBUTION TO THE RADIATED FIELD 
C 


Gall 
PLANEW(NWIRES,NW1,NW2,NP1,NP2,XH,YH,ZH,THR,P 
HR,RW) 

DO 210 L=1,MT2 

ETW=ETW+RW(L)*CNT(L+NS) 

210 EPW=EPW+RW(L+MT2)°CNT(L+NS) 
Cc 
c BOR SURFACE CONTRIBUTION TO THE RADIATED 
FIELD 
Cc 

DO 300 M=1,MHI 

NM=M- 1 

EXP1=CEXP(CMPLX(0.,FLOAT(NM)*PHR)) 

EXP2=CONJG(EXP1) 

CALL PLANES(NM,.NP1,0,NT,RH,ZH,XT,AT,THR,RS) 

NTOP 1=MODES-NM 

NTOP2=NBLOCK-(NTOP1+1) 

NS2=NTOP1°N 

NS1=NTOP2°N 

DO 250 L=1,MT1 

ET 1=ET14+RS(L)°CNT(L+NS1)*EXP1 

EP1=EP1+RS(L+N)*CNT(L+NS1)°EXP1 

IF(NM.EQ.0) GO TO 250 

ET1=ET14+RS(L)*CNT(L+NS2)"EXP2 

EP1=EP1-RS(L+N)°CNT(L+NS2)"EXP2 

250 CONTINUE 

DO 260 L=1,MP1 

ET2=ET2+RS(L+MT1)°CNT(L+NS1+MT1)°EXP1 

EP2=EP2+RS(L+MT1+N)*CNT(L+NS14+MT1)°EXP1 

IF(NM.EQ.0) GO TO 260 

ET2=ET2-RS(L+MT1)°CNT(L+NS2+MT1)°EXP2 

EP2=EP2+RS(L+MT14+N)*CNT(L+NS2+MT1)°EXP2 

260 CONTINUE 
300 CONTINUE 

ET=UC*(ET1+ET2+ETW) 

EP=UC*(EP1+EP2+EPW) 

SUM=SUM+(CABS(ET)**2+CABS(EP)**2)*STRA 

1000 CONTINUE 
2000 CONTINUE 

PRAD=T1°P1°SUM 

GAIN=4."PI°EMAX**2/PRAD 

GDB=10.*alog10(GAIN) 

WRITE(8,8010) PRAD,EMAX,GAIN,GDB 

8010 FORMAT(//,10X, TOTAL RADIATED 
POWER=’,E15.8,/,10X, 

* MAXIMUM FIELD VALUE (V/M)=',E15.8,/,10X, 

* 'GAIN=',E15.8,/,10x,'(dB)=",f8.2) 

WRITE(6,*) *PRAD,GAIN,dB=",PRAD,GAIN,GDB 

9999 STOP 

END 

SUBROUTINE 
PLANES(M1,NP1,NP2,NT,RH,ZH,XT,AT,THR,R) 

E 
C PLANE WAVE EXCITATION VECTOR FOR BOR 
ELEMENTS. 
C NO CHANGE FROM HARRINGTON'S ORIGINAL 
PROGRAM OTHER THAN 
C ** ONLY ONE MODE PER CALL ** 
€ 
COMPLEX 
R(1800),U,U1,UA,UB,FA(50),FB(50),F2A,F2B,F1A,F1B 
COMPLEX U2,U3,U4,U5,CMPLX 
DIMENSION RH(1000),ZH( 1000) ,XT(10),AT(10) 
DIMENSION R2(10),22(10),BU(5000) 
U=(0.,1.) 
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U1=3.141593°U°"M1 
NF=1 
M2=M1 
NP=NP14NP2 
MP=NP-1 
MT=MP-1 
N=MT+MP 
N2=2°N 
CC=COS(THR) 
SS=SIN(THR) 
M3=M1+1 
M4=M2+3 
IF(M1.EQ.0) M3=2 
M5=M1+2 
M6=M2+2 
DO 12 IP=1,MP 
K2=IP 
l=IP+1 
DR=.5*(RH(I)-RH(IP)) 
DZ=.5°(ZH(I)-ZH(IP)) 
D1=SQORT(DR*DR+D2Z°DZ) 
R1=.25*(RH(l)+RH(IP)) 
IF(ABS(R1).LT.1.E-5) R1=1. 
21=.5°(ZH(l)+ZH(IP)) 
DR=.5*°DR 
D2=DR/R1 
DO 13 C=1,NT 
R2(L)=R1+DR°XT(L) 
Z2(L)=Z1+DZ°XT(L) 
13 CONTINUE 
D3=DR*CC 
D4=-DZ°SS 
B5=01"°CC 
DO 23 M=M3,M4 
FA(M)=0. 
FB(M)=0. 
23 CONTINUE 
DO 15 L=1,NT 
X=ss R2(L) 
IF(X.GT..5E-7) GO TO 19 
DO 20 M=M3,M4 
BU(M)=0. 
20 CONTINUE 
Ba Zi 
S=1- 
GOTO 18 
19 M=2.8°X+14.-2./X 
IF(X.LT..5) M=11.8+ALOG10(X) 
IF(M.LT.M4) M=M4 
BJ(M)=0. 
JM=M-1 
BJU(JM)=1. 
DO 16 J=4,M 
J2=JM 
JM=JM-1 
J1=JM-1 
BU(JM)=J1/X*BU(J2)-BU(JM+2) 
16 CONTINUE 
So) 
IF(M.LE.4) GO TO 24 
DO 17 J=4,M,2 
$=S+BJ(J) 
17 CONTINUE 
24 S=BU(2)4+2.°S 
18 ARG=Z2(L)°CC 
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UA=AT(L)/S*CMPLX(COS(ARG),SIN(ARG)) 
UB=XT(L)*UA 
DO 25 M=M3,M4 
FA(M)=BJ(M)°UA+FA(M) 
FB(M)=BU(M)°UB+FB(M) 
25 CONTINUE 
15 CONTINUE 
IF(M1.NE.0) GO TO 26 
FA(1)=- FA(3) 
FB(1)=- FB(3) 
26 VA=U1 
DO 27 M=M5,M6& 
M7=M-1 
M8=M+1 
F2A=UA"(FA(M8)+FA(M7)) 
F2B=UA*(FB(M8)+FB(M7)) 
UB=U"UA 
F1A=UB*(FA(M8)-FA(M7)) 
F1B=UB°(FB(M8)-FB(M7)) 
U4=D4°UA 
U2=D3°F 1A+U4°FA(M) 
U3=D3°F 1B+U4*FB(M) 
U4=DR°F2A 
U5=DR°F2B 
K1=K2-1 
K4=K14+N 
K5=K2+N 
R(K2+MT)=-D5*(F2A+D2°F2B) 
R(K5+MT)=D1°(F 1A+D2°F 1B) 
IF(IP.EQ.1) GO TO 21 
R(K1)=R(K1)+U2-U3 
R(K4)=R(K4)4U4-U5 
IF(IP.EQ.MP) GO TO 22 
21 R(K2)=U2+U3 
R(K5)=U44+U5 
22 K2=K2+N2 
UA=UB 
27 CONTINUE 
12 CONTINUE 
RETURN 
END 
SUBROUTINE 
PLANEW(NWIRES,NW1,NW2,NP1,NP2,XH,YH,ZH, 
* THR,PHR,R) 
C 
C PLANE WAVE EXCITATION VECTOR ELEMENTS FOR 
WIRE AND 
C ATTACHMENT SEGMENT. INCIDENCE DIRECTION IS 
(THR,PHR). 
C -ONLY ONE INCIDENCE ANGLE PER CALL 
C -WIRE DOES NOT HAVE TO LIE IN THE Z=0 PLANE 
C 
COMPLEX 
U0,AA,BB,C,R(400),CEXP,EXP,F11,Fl2,S!,.DI,CMPLX 
DIMENSION 
ZH(400),NS(2),NW1(2),NW2(2),XH(400),YH(400) 
MP2=NP2-1 
MT2=NP2-2 
MT22=2°MT2 
DO 5 L=1,NWIRES 
5 NS(LJ=ENW2(L)-NW1(L)+1 
U0=(0.,0.) 
CC=COS(THR) 
SS=SIN(THR) 
CP=COS(PHR) 


SP=SIN(PHR) 
UP=Ss cr 
VP=SS*SP 
WP=CC 
DO 12 iP=1,MP2 
t=tP+NP1 
l=tt+1 
ZS=.5*(ZH(!)+ZH(It)) 
XS=.5*(XH(l)+XH(Ii)) 
YS=.5*(VH(i)+YH(Il)) 
DX=XH(I)-XH(th) 
DY=YH(I)- YH (It) 
D1=SQRT(DX**2+DY**2) 
SU=DY/D1 
CU=DX/D1 
C FOR WIRES IN PERPENDICULAR TO Z SIN(V)=1 AND 
COS(V)=0 
SV=IL0 
CV=0.0 
C 
C WIRE SEGMENT CALCULATIONS 
C 
A=UP*CU+VP*SU 
B=UP*XS+VP*YS+WP°ZS 
C=CMPLX(0.,A) 
EXP=CEXP(CMPLX(0.,B)) 
AA=CC*(CU°CP+SU*SP) 
BB=SU*CP-SP*CU 
PS!l=D1°A/2. 
IF(PSt.NE.0.) GO TO 60 
SINC=s1. 
GOTO 61 
60 SINC=SIN(PSI)/PSI 
61 COSP=COS(PSI) 
FI1=SINC*D1°EXP/2. 
F\2=(0.,0.) 
IF(ABS(A).LT.1.E-4) GO TO 62 
CSP=COsr-siNe 
IF(ABS(CSP).LT.1.E-4) GO TO 62 
FI2=EXP/C°CSP 
62 CONTINUE 
SI=FI1+Fi2 
DI=Ft1-Fi2 
C 
C R-WIRE-THETA 
C 
lIF(tP.EQ.MP2) GO TO 10 
R(IP)=AA*S! 
R(tP+MT2)=BB°SI 
10 CONTINUE 
Cc 
C R-WIRE-PHI 
C 
14 IF(IP.EQ.1) GO TO 12 
R(IP-1)=R(tP-1)+AA*D! 
R(iP-1+4MT2)=R(IP-1+MT2)+BB°Di 
12 CONTINUE 
RETURN 
END 


Reflector 
C PROGRAM mrc.f (Marine Radio Antenna Program) 


C 
¢ Modified from ORDONEZ.f -- hyperboloid becomes a disc 


C Paraboloie WITH A PARASITIC RING DIPOLE FEED. 
Fed dipole is at z=0. 
C This program uses a parasitic dipole and two h-plane 
C parasitic dipoles and the subreflector. This system also 
contains 
C acavity ring. 
C 
C ICALC=0 CALC CURRENTS AND FIELD 
Cc =1 CALC FIELDS ONLY (READ CURRENTS) 
C IMP=0 PEREECTLY CONDUCTING SURFACES 
Cc =1 SOME NONZERO SURFACE RESISTANCE 
C ICWRT=0 WRITE CURRENT COEFFICIENCTS TO 
FILE pceurrent 
C tRES=0 READ RESISTIVE CORRECTIONS 
Cc 
COMPLEX 
EP,ET,Z(500000),RS(1000),RVW(400),B(700),C(700),U UC 
COMPLEX ZL0(400),ZL(1500),EC,EX,ZSS(50000) 
COMPLEX EXP1,EXP2,ZSW(50000), ZWW(10000) 
COMPLEX 
CEXP,CONJG,CMPLX,ET1,ET2,EP1,EP2, ETW,EPW 
DIMENSION 
RH(400),ZH(400),XT(4),AT(4),IPS(700), NW1(4),NW2(4) 
DIMENSION 
A(400),X(400),EXP(800), ANG(800),ECP(800) 
DIMENSION ECV(800),EXV(800),PHC(800),PHX(800) 
DIMENSION 
XX(41), YY¥(41),XH(400), YH(400),PHA(400) 
DATA 
PILSTART ,STOP,DT/3.14159,0.,180.,1./, \CWRT/0/,IRES/1/ 
DATA ICALC/0/,IPRINT/O/,IMP/0/, ITEST/1/ 
DATA nt/2/,xt(1),at(1)/.5773503,1./ 
¢ data nt/4/,xt(1),xt(2),at(1),at(2)/.33998104,.8611363115, 
* 6521451548, .3478548451/ 
OPEN(1,FILE='outgaus’,status='old’) 
OPEN(30,FiLE='rh.m’) 
OPEN(31,FILE='zh.m’) 
READ(1,°) NPHI 
DO 3 K=1,NPHI 
READ(1,*) X(K),A(K) 
3 CONTINUE 
RAD=P!/180. 
BK=2.*Pl 
U=(G7 1) 
U0=(0.,0.) 
UC=-U/4./Pi 
NT2=NT/2 
DO 1 K=1,NT2 
K1=NT-K+1 
AT(K1)=AT(K) 
XT(K1)=XT(K) 
1 XT(K)=-XT(K) 
C 
C NP10 = NUMBER OF PARABOLOID POINTS 
C NP11 =NUMBER OF DISC POINTS 
C NP12 = NUMBER OF RING POINTS 
C 
MODES=1 
MHI=MODES+1 
NBLOCK=2*MODES +1 
NWIRES=1 
WRITE(6,*) "ENTER START,STOP,DT’ 
READ(5,*) START,STOP,DT 
NP10=61 
NP11=31 
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NP12=3 
WRITE(6,°) ‘ENTER PARA (NP 10), DISC (NP11) AND 
RING (NP12) PTS' 
READ(5,*) NP10,NP11,NP12 
NP1=NP104NP11+NP12 
Cc 
C DM=MAIN REFLECTOR DIAMETER 
C DS=DISC (SUBREFLECTOR) DIAMETER 
C ZC=DISC LOCATION ALONG Z AXIS 
C FOD=f/OM RATIO 
C RC=RADIUS OF THE RING 
C SW=RING (STRIP) WIDTH 
Cc 
DM=5.0 
RO=DM/2. 
DS=1.1 
ZC=.25 
FOD=.3 
RC=0.5 
SWe.05 
WRITE(6,°) ‘ENTER DM, DS, ZC (>0), FOD, RC, ZR 
AND SW' 
READ(5,*) OM,OS,ZC,FOD,RC,ZR,SW 
FM=FOD*DM 
WRITE(6,°) ‘ENTER ZSHIFT FOR PARABOLOID' 
READ(5,°) ZSHIFT 
C 
C GENERATE PARABOLOID CONTOUR 
Cc 
PHIV=2.°ATAN(1./4./FOD) 
FNM=FLOAT(NP 10-1) 
DO 50 l=1,NP10 
TH=FLOAT(I-1)*PHIV/FNM 
RM=2.°FM/(1.4COS(TH)) 
XH(I)=0. 
TAO)=0. 
ZH(l)=-RM*COS(TH)*BK+ZSHIFT “BK 
RH(l)=RM°SIN(TH)*BK 
PHA(II)=0. 
50 CONTINUE 
Cc 
C DISC LOCATED A Z=ZC 
Cc 
OSK=DS/2./FLOAT(NP11-1) 
DO 5 Il=1,NP11 
N=l4NP10 
XH(Il)=0. 
YH(I)=0. 
ZH(I)=ZC°BK 
RH(iIJ=FLOAT(I-1)*DSK*BK 
PHA(II)=0. 
5 CONTINUE 
Cc 
C NEXT THE SIDEWALLS 
Cc 
DSW=SW/FLOAT(NP 12-1) 
ZRBK=ZR°*BK 
IF(NP12.NE.0) THEN 
DO 7 l=1,NP12 
Il=NP11+NP10+I 
XH(Il)=0. 
YH(IN=0. 
RH(Il)=RC*BK 
ZH (Il)=FLOAT(I-1)°DSW*BK+ZRBK 
7 PHA(I=0. 
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ENDIF 
NP1=NP104NP114NP12 
C 
C DIPOLE ANTENNA POINTS IN WAVELENGTHS: AW= 
RADIUS; NP2= NUMBER 
C OF POINTS (=NP21+NP22; NP21 MUST BE ODD); 
DLGTH= LENGTH; 
C PLGTH IS PARASITIC ELEMENT LENGTH. 
C 
AWz=.001 
DLGTH=.45 
PLGTH=.44 
WRITE(6,”) ‘ENTER DLN, PLN, DDIP, ROT AND AW’ 
WRITE(6,*) '(DDIP<0 IS TOWARD PARABOLOID)’ 
READ(5,*) DLGTH,PLGTH,DDIP,ROT,AW 
CENT=DLGTH/2. 
AK=AW'BK 
NP21=11 
NP22=0 
WRITE(6,°) ‘ENTER NP21 AND NP22 (ONLY 1 DIP IF 
NP22=0)’ 
READ(5,*) NP21,NP22 
NW1(1)=NP 141 
NW2(1)=NP1+NP21 
NP2=NP214+NP22 
DEL=DLGTH/FLOAT(NP 21-1) 
DO 10 l=1,NP21 
ll=NP1+1 
YRHUD=0. 
ZH (I)=0. 
XH(Ih=BK*(FLOAT(I-1)°DEL-CENT) 
RH(I)=SQRT (XH (II) °°24+YH(I1)**2) 
PHA(IN=ATAN2(YH(II),XH(I1)+1.E-6) 
10 CONTINUE 
IF(NP22.NE.0) THEN 
NWIRES=2 
NW1(2)=NP14+NP21+1 
NW2(2)=NP 14NP214NP22 
DEL=PLGTH/FLOAT(NP22-1) 
CENT=PLGTH/2. 
BO 1W1W=1,.NP2Z2 
Il=NP14+NP214! 
ZH(II}=DDIP*BK 
XH(Il)=BK*(FLOAT(I-1)*DEL-CENT)*COS(ROT*RAD) 
YH(I)=BK*(FLOAT(I-1)°DEL-CENT)*SIN(ROT*RAD) 
RH(IN)=SQRT (XH(I1)°°24+Y HII) °*2) 
PHA(IN=ATAN2(YH (II) XH(11)+1.E-6) 
11. CONTINUE 
ENDIF 
C NP23=NUMBER OF H-PLANE DIPOLE POINTS 
C HLGTH=H-PLANE DIPOLE LENGTH 
C DHP=H-PLANE DIPOLE DISTANCE FROM FEED 
C ZHP=H-PLANE DIPOLE Z COORDINATE 
C NP24=NUMBER OF POINTS FOR THE OTHER 
H-PLANE DIPOLE 
NP23=7 
NP24=7 
HLGTH=.2 
DHP=.2 
ZHP=.001 
WRITE(6,°) ‘ENTER NUMBER OF POINTS PER 
H-PLANE DIPOLE’ 
WRITE(6,*) ‘(0 FOR NO H-PLANE PARASITIC 
DIPOLES)’ 
READ(5,*) NP23 


IF(NP23.NE.0) THEN 

WRITE(6,°) 'H-PLANE DIPOLE LENGTH, DISTANCE 
AND LOCATION’ 

READ(5,*) HLGTH,DHP,.ZHP 

NWIRES=NWIRES+2 

NP24=NP23 

NP2=NP214+NP224+NP234+NP24 

NW1(3)=NP14+NP214+NP22+1 

NW2(3)=NP14NP214+NP22+NP23 

NW1(4)=NP14+NP214NP22+NP2341 

NW2(4)=NP 14NP214NP224+NP234+NP24 

DEL=HLGTH/FLOAT(NP23-1) 

CENT=HLGTH/2. 

DOWZ I= Ne2zs 

l11=NP14+NP214+NP22+1 

ZH(I1)=ZHP*BK 

YH(l1)=DHP*BK 

XH(I1)=BK*(FLOAT(I-1)°DEL-CENT) 

RH(11)=SQRT(XH(11)°°2+YH(11)°*2) 

PHA(I1)=ATAN2(YH(I1),XH(11)+1.E-6) 

I2=11+4NP23 

ZH(I2)=ZHP*BK 

YH(I2)=-DHP*BK 

XH(I2)=BK*(FLOAT(I-1)*DEL-CENT) 

RH(I2)=SQRT(XH(12)**2+YH(I2)**2) 

PHA(I2)=ATAN2(YH(I2),XH(I2)+1.-6) 

12 CONTINUE 

ENDIF 

CLOSE(1) 

NP=NP1+NP2 

WRITE(6,*) "ENTER FREQUENCY FACTOR’ 

READ(5,*) FFAC 

OPEN(8 file=’outmrc’) 

WRITE(8,8000) 
DM,DS,FOD,PHIV/RAD,NP10,NP11,RC,SW, 
NP12,ZSHIFT,OLGTH,PLGTH,DDIP,NP2,AW,MODES,NT, 
NPHI,FFAC 

IF(NP23.NE.0) WRITE(8,8001) HLGTH,DHP,ZHP,NP23 

8000 FORMAT(//,5X,"*** MARINE RADIO REFLECTOR 
ANTENNA °°**' //,5xX, 

“REFLECTOR PARAMETERS (LENGTHS IN 
WAVELENGTHS):’//,5X, 

*'MAIN REFL DIAz’ F8.3,/,5X,'DISC 
DIA=',F8.3,/,5X,'F/D=',F8.4, 

*/,5X,"PHIV (DEG)=',F8.2,/,5X,,MAIN REFL POINTS 
(NP11)=",14,/, 

* 5X,,SUB REFL POINTS (NP12)=',14,/,5X,/,5X, 

* ‘RING RADIUS=’,F8.3,/,5X,’SIDEWALL 
LENGTH=’,F8.3,/,5X, 

* 'NUMBER OF RING POINTS (NP12)=',14,//,5X, 

**PARABOLOID SHIFT DISTANCE, 
ZSHIFT=',F8.3,/,5X, 

*'DIPOLE PARAMETERS:’,/,5X,'FED DIPOLE 
LENGTH='F8.3,/,5X, 

* PARASITIC DIPOLE 
LENGTH=',F8.3,/,5X,,SPACING=',F8.3,/,5X, 

*'NUMBER OF DIPOLE POINTS =’',14,/,5X,’RADIUS OF 
DIPOLE=’, 

*F8.4,//,5X, NUMBER OF AZIMUTHAL 
MODES=',13,/,5X,'NT=',13, 

* /,5X,'NPHIl=',13,/,5X,,F REQ. SCALE FACTOR=’ F8.4) 

8001 FORMAT(/,5X,’H-PLANE DIPOLE 
LENGTH=',F8.3,/,5X, 


* 'H-PLANE DIPOLE DISTANCE FROM 
AXIS=',F8.3,/,5X, 

* 'Z POSITION OF H-FIELD — 
DIPOLE=',F8.3,/,5X,,NUMBER OF H-FIELD 

* DIPOLE POINTS="',13) 

IF(ISEG.EQ.0) WRITE(8,1300) 

1300 
FORMAT(/,5X,"INDEX’,8X,'Z(I)',10X,’RHO(1)',8x,'ZSURF’) 
c 
C FILE difdat CONTAINS RESISTIVE CORRECTIONS IF 
IRES.NE.O 
c 

IF((IMP.EQ.0).OR.(IRES.NE.0)) GO TO 70 
OPEN(9,FILE='difdat’) 
READ(Q,") ID 
WRITE(6,*) ‘ID=',ID 
DO 75 l=1,ID 
READ(9,*) II,RI,YY(1),DUM 
IF(YY(I).LT.O.) YY(I)=0. 
75 CONTINUE 
XXl=1./FLOAT(ID-1) 
DO 15 Il=1,ID 
XX(Il)=XXIFLOAT(II-1) 
15 CONTINUE 
WRITE(6,*) “SURFACE IMPEDANCE VALUES READ' 
70 CONTINUE 
DO 52 I=1,NP1 
IF(ABS(ZH(I)).LT..001) ZH(I)=0. 
IF(ABS(RH(I)).LT..001) RH(1)=0. 
ZL0(1)=(0.,0.) 
XH(I)=XH(I)*FFAC 
YH(I)=YH(l)*FFAC 
ZH(I)=ZH(1)"FFAC 
RH(I)=RH(I)"FFAC 
ZHB=2H(I)/BK 
RHB=RH(1)/BK 
WRITE(30,5019) RHB 
WRITE(31,5019) ZHB 
¢  f(rhb.ge.10.) zIO(i)=(.2,0.) 
Ec 
C NONZERO SURFACE IMPEDANCES IF IMP > 0 
Cc 
IF((IMP.EQ.0).OR.(IRES.NE.0)) GO TO 51 
IF(I.GT.NP10) GO TO 51 
Cc 
C USE THE EXACT RESISTIVE CORRECTIONS -- 
INTERPOLATE BETWEEN 
C DATA POINTS. 
e 
UU=RHB/RO 
CALL INTERP(UU,VV,XX, YY, ID) 
ZLO(I)=VV 
51 CONTINUE 
WRITE(8,8004) I,ZHB,RHB,ZL0(1) 
52 CONTINUE 

8004 FORMAT(6X.14,4X,F8.3,8X,F8.3,6X ,2F8.4) 
WRITE(8,1310) 

1310 FORMAT(/,5X,'DIPOLE COORDINATES: ’,//, 
* 5X," INDEX’ ,8X,'Z()’,10X,’RHO(I)',10X,’PHI PLN’) 
DO 53 I=NP14+1,NP 
IF(ABS(ZH(I)).LT..001) ZH(I)=0. 
IF(ABS(RH(1)).LT..001) RH(1)=0. 
XH(I)=XH (I) *FFAC 
YH(I)=YH(I)"FFAC 
ZH(I)=ZH(l)*FFAC 
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RH(!)=RH(l)*FFAC 

PHB=PHA(1I)/RAD 

ZHB=ZH(1)/BK 

RHB=RH(1)/BK 

WRITE(30,5019) RHB 

WRITE(31,5019) ZHB 

53 WRITE(8,8005) I|,ZHB,RHB,PHB 
8005 FORMAT(6X,14,4X,F8.3,8X,F8.3,9X,F8.2) 
799 CONTINUE 

IF(ITEST.EQ.0) GO TO 9998 

write(6,") ‘geometry defined’ 

Cc 

C DEFINE DIMENSIONS OF THE IMPEDANCE MATRIX 

BLOCKS 

Cc 
MT1=NP1-2 
MP1=NP 1-1 
N=MP1+MT1 
MT2=NP2-2 
NSURF=NBLOCK*'N 
NSW=NSURF+MT2 
NROW=NSW 
WRITE(6,°) 

‘NP 1,NP2,NSURF,NROW=',NP1,NP2,NSURF,NROW 
WRITE(6,*) 'MT2 SHOULD BE ODD -- MT2=',MT2 
IF(ITEST.EQ.0) GO TO 9998 
IF(ICALC.EQ.0) THEN 

Cc 

C COMPUTE IMPEDANCE MATRIX ELEMENTS 

Cc 
CALL 

ZMATWW(NWIRES,NW1,NW2,XH,YH,RH,ZH,NT,.XT,AT.A 

K,ZWW) 

CALL ZASMBO(NP1,NP2,MODES,Z,ZWW) 

write(6,°) ‘wire impedance computed’ 

IF(IMP.NE.0) THEN 

CALL ZLOAD(NP1,RH,ZH,ZL0,ZL) 

WRITE(6,*) ‘RESISTIVE SURFACE IMPEDANCE 
MATRIX COMPUTED’ 

ENDIF 

DO 400 M=1,MHI 

NM=M- 1 

CALL 
ZMATSS(NM,NM,NP10,NP11,NP12,NPHI,NT,RH,ZH,X,A,X 
T,AT,ZSS) 

IF(IMP.NE.O) CALL ZTOT(MT1,MP1,ZL.ZSS) 

CALL 
ZMATSW(NWIRES,NW1,NW2,NP10,NP11,NP12,XH,YH,R 
H,ZH, 

* NT .XT,AT,NPHI,X,A,NM,AK,ZSW) 

CALL ZASMBN(NP1,NP2,MODES,NM,Z,Z2SS.ZSW) 

IF(NM.EQ@.0) GO TO 400 

NMN=-NM 

CALL 
ZMATSW(NWIRES NW1,NW2,NP10,NP11,NP12,XH,YH,R 
FZ 

* NT .XT,AT,NPHI.X,A,NMN,AK,ZSW) 

CALL ZASMBN(NP1,NP2,MODES,NMN.Z,ZSS.Z2SW) 

400 CONTINUE 

write(6,°) 'z filled’ 

CALL DECOMP(NROW,IPS.Z) 

write(6,*) 'z decomposed’ 

Cc 

C EXCITATION ELEMENTS (NONZERO FOR THE 

DIPOLE ONLY) 


Cc 
DO 40 l=1,NROW 
40 B(I)=(0.,0.) 
JO=NSURF+(NP21-1)/2 
B(!0)=(1.,0.) 
CALL SOLVE(NROW, IPS,Z,B,C) 
write(6,*) ‘current coefficients computed’ 
c do 665 i=1,nrow 
c665 write(6,°*) ii,c(ii) 
IF((ICWRT.EQ.0). AND.(ICALC.EQ.0)) THEN 
Cc 
C WRITE CURRENTS ON DISC FOR PATTERN 
INTEGRATION 
c (program ‘cavarchint.f’) 
Cc 
OPEN(S file='pcurrent’) 
WRITE(3,°) 
DM,AW,ZC,RC,DS,NP1,NP2,NBLOCK,NWIRES ,FFAC 
WRITE(3,°) 
NP,MT1,MP1,N,MT2,NSURF NROW,MHI,MODES 
WRITE(3,*) (XH(1),l=1,NP) 
WRITE(3,*) (YH(I),l=1,NP) 
WRITE(S3,*) (ZH(I),l=1,NP) 
WRITE(3,*) (RH(I),l=1,NP) 
WRITE(3,*) (PHA(I),l=1,NP) 
WRITE(3,*) (NW1(1),l=1,.NWIRES) 
WRITE(3,*) (NW2(1),l=1,NWIRES) 
WRITE(3,*) (C(I) ,l=1,NROW) 
ENDIF 
ENDIF 
Cc 
C iF ICALC.NE.0 THEN READ CURRENT COEFFICIENTS 
FROM DISK FILE 
Cc 
IF(ICALC.NE.0) THEN 
OPEN(3 file='pcurrent’*, STATUS="OLD’) 
READ(3,*) 
DMX, AWX,ZCX,RCX.DSX,NP1X,NP2X ,NBLX,NWX 
READ(3,*) 
NPX,MT 1X ,MP1X,NX,MT2X ,NSX,NRX,MHIX MODESX 
READ(3,*) (xh(i),i=1,np) 
READ(3,*) (yh(i),i=1,np) 
READ(3,*) (zh(i),i=1,np) 
READ(3,*) (rh(i),i=1,np) 
READ(3,*) (pha(i),i=1,np) 
READ(3,*) (nw1(i),i=1,nwires) 
READ(3,*) (nw2(i),i=1,nwires) 
READ(3,*) (c(i) ,i=1,nrow) 


CLOSE(3) 
write(6,*) 'data read from pcurrent' 
do 909 i=1,np 

909 write(6,*) ‘i,rh,zh=',i,rh(i)/ok,zh(i)/ok 
ENDIF 
IT=INT((STOP-START)/DT)+1 


Cc 
C RECEIVER PHI CUTS: DO 0 AND 90 (E- AND H- 
PLANES) 
C 
DO 501 IP=1,2 
ECX=1.e-10 
PHRO=0. 
IF(IP.EQ.2) PHRO=RAD‘90. 
DO 500 I=1,IT 
THETA=FLOAT(I-1)*DT+START 
THR=THETA‘RAD 
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PHR=PHRO 
IF(THETA.LE.180.) GO TO 99 
THR=(360.-THETA)*RAD 
PHR=PHRO+PI 
99 CONTINUE 
c  write(6,°) ‘theta, phi=’ theta,phr/rad 
ET1=U0 
ET2=U0 
EP1=U0 
EP2=U0 
ETW=U0 
EPW=U0 
C 
C DIPOLE FIELD 
Cc 
CALL 
PLANEW(NWIRES,NW1,NW2,NP1,NP2,XH,YH,ZH,THR,P 
HR,RW) 
DO 210 L=1,MT2 
ETW=ETW-+RW(L)°C(L+NSURF) 
210 EPW=EPW+RW(L+MT2)*C(L+NSURF) 
C 
C REFLECTOR AND CAVITY FIELD CONTRIBUTIONS 
Cc 
DO 300 M=1,MHI 
NM=M-1 
EXP1=CEXP(CMPLX(0.,FLOAT(NM)*PHR)) 
EXP2=CONJG(EXP1) 
CALL 
PLANES(NM,NP10,NP11,NP12,NT,RH,ZH,XT,AT,THR,RS) 
NTOP1=MODES-NM 
NTOP2=NBLOCK-(NTOP1+1) 
NS2=NTOP1‘N 
NS1=NTOP2'‘N 
DO 250 L=1,MT1 
ET1=ET1+RS(L)°C(L+NS1)"EXP1 
EP1=EP1+RS(L+N)*C(L+NS1)*EXP1 
IF(NM.EQ.0) GO TO 250 
ET 1=ET1+RS(L)*C(L+NS2)*EXP2 
EP1=EP1-RS(L+N)*C(L+NS2)*EXP2 
250 CONTINUE 
DO 260 L=1,MP1 
ET2=ET2+RS(L+MT1)°C(L+NS14+MT1)*EXP1 
EP2=EP2+RS(L+MT 14+N)*C(L+NS14+MT1)"EXP1 
IF(NM.EQ.0) GO TO 260 
ET2=ET2-RS(L+MT 1)*C(L+NS24+MT1)*EXP2 
EP2=EP2+RS(L+MT14+N)°C(L+NS2+MT1)*EXP2 
260 CONTINUE 
300 CONTINUE 
ET=UC*(ET1+ET2+ETW) 
EP=UC*(EP1+EP2+EPW) 
EC=ET 
EX=EP 
ECV (I)=CABS(EC) 
EXV(Il)=CABS(EX) 
ECR=REAL(EC) 
ECI=AIMAG(EC) 
EXR=REAL(EX) 
EXl=AIMAG(EX) 
PHC(I)=ATAN2(ECI,ECR+1.e-10)/RAD 
PHX(I)=ATAN2(EXI,EXR+1.e-10)/RAD 
ANG(I)=THETA 
ECX=AMAX1(ECX,ECV(l),EXV (I) 
500 CONTINUE 
WRITE(8,103) PHR/RAD,ECX 


103 FORMAT(/,5X, PHI OF RECEIVER 
(DEG)=’,F8.2,/,5Xx, 
* "MAXIMUM FIELD VALUE (V/M)=',E15.5) 
DO 600 l=1,IT 
ECP(I)=(ECV(I)/ECX)**2 
EXP(I)=(EXV(I)/ECX)**2 
ECP(I)=AMAX1(ECP(I),.00001) 
EXP(I)=AMAX1(EXP(1!),.00001) 
ECP(l)=10.°ALOG10(ECP(I)) 
EXP(l)=10.°ALOG10(EXP(I)) 
600 CONTINUE 
IF(IPRINT.NE.O) GO TO 310 
WRITE(8,5015) 
5015 
FORMAT (///,7X,,ANGLE’,15X,’CO-POL’,.25X,K-POL',/,7X, 
*(DEG)',4X,2((VOLTS)’,4X,'(DEG)',3X,’(DB-REL)',4X)) 
DO 9000 L=1,IT 
WRITE(8,5016) 
ANG(L),ECV(L),PHC(L),ECP(L),EXV(L),PHX(L) 
* EXP(L) 
5016 FORMAT(5X,F6.1,3X ,2(F8.4,3X,F7.1,3X,F7.2,3X)) 
9000 CONTINUE 
310 CONTINUE 
Cc 
C WRITE MATLAB FILES ( .m FILES) 
Cc 
IF(IP.EQ.1) THEN 
OPEN(2,file='ang.m’) 
OPEN(3 ,file=’cpole.m’) 
OPEN(4 ,file=’xpole.m’) 
DO 9097 l=1,IT 
WRITE(2,5019) ANG(I) 
WRITE(3,5019) ECP(1) 
9097 WRITE(4,5019) EXP(l) 
CLOSE(2) 
CLOSE(3) 
CLOSE(4) 
ENDIF 
IF(IP.EQ.2) THEN 
OPEN(3,file=’cpolh.m’) 
OPEN(4 file=xpolh.m’) 
DO 9098 [=1,IT 
WRITE(3,5019) ECP(I) 
9098 WRITE(4,5019) EXP(I) 
CLOSE(3) 
CLOSE(4) 
ENDIF 
501 CONTINUE 
5019 FORMAT(F8.3) 
9998 STOP 
END 
SUBROUTINE SOLVE(N,IPS,UL,B,X) 
COMPLEX UL(500000),B(700),X(700),SUM 
DIMENSION IPS(700) 
NP 1=N+1 
IP=IP st) 
X(1)=B(IP) 
DO 2 l=2,N 
IP=IPS(1) 
IPB=IP 
IM1=I-1 
SUM=(0.,0.) 
DO 1 J=1,IM1 
SUM=SUM+UL(IP)*X(J) 
1 IP=IP+N 


2 X(I)=B(IPB)-SUM 
K2=N°(N-1) 
IP=IPS(N)+K2 
X(N)=X(N)/UL(IP) 
DO 4 IBACK=2,N 
lI=NP1-IBACK 
K2=K2-N 
IPI=IPS(1)+K2 
IP 1=1+1 
SUM=(0.,0.) 
IP=IPI 
DO 3 J=IP1,N 
IP=IP+N 
3 SUM=SUM+UL(IP)*X(J) 
4 X(N=(X(1)-SUM)/UL(IPD 
RETURN 
END 
SUBROUTINE DECOMP(N,IPS,UL) 
COMPLEX UL(500000) ,P!VOT,EM 
DIMENSION SCL(700),IPS(700) 
DO 5 Il=1,N 
IPS(I=! 
RN=0. 
=| 
OS 2J=1,N 
ULM=ABS(REAL(UL(J1)))+ABS(AIMAG(UL(J1))) 
J1=J14+N 
IF(RN-ULM) 1,2,2 
1 RN=ULM 
2 CONTINUE 
SCL(l=1./RN 
5 CONTINUE 
NM1=N-1 
K2=0 
DO 17 K=1,NM1 
BIG=0. 
DO 11 l=K,N 
{\P=IPS(l) 
IPK=!IP+K2 


SIZE=(ABS(REAL(UL(IPK)))+ABS(AIMAG(UL(IPK))))°"SCL( 
IP) 

IF(SIZE-BIG) 11,11,10 

10 BIG=SIZE 
IPV=! 

11 CONTINUE 
IF(IPV-K) 14,15,14 

14 J=!IPS(K) 
IPS(K)=IPS(IPV) 
IPS(IPV)=J 

15 KPP=IPS(K)+K2 
PIVOT=UL(KPP) 
KP 1=K+1 
DO 16 I=KP1,N 
KP=KPP 
IP=IPS(I)+K2 
EM=-UL(IP)/PIVOT 

18 UL(IP)=-EM 
DO 16 J=KP1,N 
IP=IP+N 
KP=KP+N 
UL(UIP)=UL(IP)+EM°UL(KP) 

16 CONTINUE 
K2=K2+N 

17 CONTINUE 


Cc 


RETURN 

END 

FUNCTION BLOG(x) 

IF(X.GT..1) GO TO 1 

X2=X°X 
BLOG=((.075*X2-.1666667)"X2+1.)°X 
RETURN 


1 BLOG=ALOG(X+SQRT(1.+X*X)) 


RETURN 
END 
SUBROUTINE ZLOAD(NP,RH,ZH,Z0,2Z) 


C COMPUTES IMPEDANCE MATRIX ELEMENTS FOR 
LOADED BODIES OF REV 

C Z0(I) IS THE SURF IMPEDANCE OF THE ITH 
SEGMENT (NP-1 SEGMENTS) 

C 2(.) ARE THE IMPEDANCE MATRIX TERMS 
(TRIDIAGONAL FOR T-T 

C SUBMATRIX; DIAGONAL FOR P-P SUBMATRIX). 
STORED IN COL VECTOR. 


Cc 


COMPLEX 


C1,C2,20(400),Z(1500),.X1,X2,X3,Y1,¥2,¥3,FN(400) 
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COMPLEX U1,U2,U3.XI,Y! 
DIMENSION RH(400),ZH(400),RS(400),0(400), SV(400) 
Pl=3.14159 

MT=NP-2 

MP=NP-1 

N=MT+MP 

DO 10 IP=2,NP 

ll=IP-1 

DR=RH(!P)-RH(I!) 
DZ=ZH(!P)-ZH(!1) 
O(Wj=SQRT(DOR*DR+D2Z°DZ) 
SV ({)=DR/D(I) 
RS(i=.5°(RHUP)+RHCII)) 
DS=D(Il)’SVi(II)/2. 


Q1=RS(1I)}+DS 
Q2=RS(I!)-DS 
FN(II)=1. 


IF((ABS(Q2).GT.1.E-6).AND.(ABS(Q1).GT.1.E-6)) 
* FN(II)=ALOG(Q1/Q2) 


10 CONTINUE 


LO=MT*3-2 

DO 20 !=1,MP 

C1=P1°Z0(!) 

IF(1.EQ.MP) GO TO 80 

Kl=2 

IF (1.EQ.1) Kl=1 

IF(ILEQ.MT) Ki=3 

=l+1 

C2=PI*Z0(II) 

A=SV(I) 

IF(ABS(A).LT.1.E-6) GO TO 41 
X1=C1°FN(1)/2./A 
X2=C1°2./A*(1.-RS(1)"FN(1)/D(1)/A) 
X3=-X2°RS(1)/D(I)/A 

GO TO 42 


41 CONTINUE 


X1=C1/2./RS(1)"D(1!) 
X2=(0.,0.) 
X3=C 1°D(1)/6./RS(!) 


42 CONTINUE 


A=SViI!) 
IF(ABS(A).LT.1.E-6) GO TO 45 


Y1=C2°FN(II)/2./A 
Y2=C2°2./A*(1.-RS(I)°FN(N)/D(II)/A) 
Y3=-Y2°*RS(II)/D(II)/A 
GO TO 40 
45 CONTINUE 
Y1=C2/2./RS(Il)*D(II) 
Y2=(0.,0.) 
Y3=C2°D(I1)/6./RS(II) 
40 CONTINUE 
Cc 
C DEFINE TRIDIAGONAL ELEMENTS FOR T-T 
SUBMATRIX (STORED IN COLS) 
C (U1- DIAG; U2- LOWER; U3- UPPER) 
Cc 
Xl=X14X24+X3 
Yl=Y1-Y24Y3 
IF(KI.EQ.1) Xl=C1/SV(1) 
C IF(KILEQ.3) Yl=C2/SV(II) 
U1=XI+YI 
U2=X 1-X3 
U3=Y1-Y3 
L=2+(I-2)°3 
IF(KI.EQ.1) L=0 
L1=L+1 
L2=L+2 
L3=L+3 
go to (50,60,70),ki 
50 Z(L1)=U1 
7& (les 0 
GO TO 80 
60 Z(L1)=U3 
2(L2)=U1 
Z(L3)=U2 
GO TO 80 
70 Z(L1)=U3 
Z(L2)=U1 
80 Z(LO+I)=2.°C 1°D(I)/RS(I) 
20 CONTINUE 
RETURN 
END 
SUBROUTINE ZTOT(MT,MP,ZL.Z) 
Cc 
C ADDS THE SURF IMPEDANCE TERMS TO THE 
TRIDIAGONAL ELEMENTS OF 
C THE BOR IMPEDANCE MATRIX Z. 
Cc 
COMPLEX 2ZL(700),2(50000) 
N=MT+MP 
MO=MT*3-2 
DO 100 l=1,MP 
LO=MT‘*N+(I-1)°N+MT 
IF(l.EQ.MP) GO TO 80 
Kl=2 
IF(1.EQ.1) Kl=1 
IF(l.EQ.MT) Kl=3 
L2=(I-1)°N+l 
L1=L2-1 
L3=L2+1 
M=2+3°*(I-2) 
IF(KI.EQ.1) M=0 
M1=M+1 
M2=M+2 
M3=M+3 
go to (50,60,70),ki 
50 Z(L2)=Z(L2)+ZL(M1) 
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2(L3)=Z(L3)+ZL(M2) 
GO TO 80 
60 Z2(L1)=2(L1)+ZL(M1) 
2(L2)=Z(L2)+ZL (M2) 
Z(L3)=Z2(L3)+ZL(M3) 
GO TO 80 
70 Z(L1)=Z(L1)+ZL(M1) 
2(L2)=2(L2)+ZL(M2) 
80 Z(L0+!)=Z2(LO+1)+ZL(MO+!) 
100 CONTINUE 
RETURN 
END 
subroutine interp(u,v,x,y,nn) 
¢ program to interpolate linearly between arrays of x and y 
c values (with dimensions nn). u is specified and v is 
returned 
¢ uniform sampling in x (the ordinate) is assumed. 
dimension x(nn),y(nn) 
delx=x(2)-x(1) 
n1=int(u/delx)+1 
n2=n1+1 
dety=y(n2)-y(n1) 
sgn=sign(1.,dety) 
dely=abs(dely) 
alpha=u-x(n1) 
zeta=dely °alpha/detx 
v=y(n1)+sgn*zeta 
return 
end 
SUBROUTINE 
PLANEW(NWIRES,NW1,NW2,NP1,NP2,XH,YH,ZH, 
* THR,PHR,R) 
Cc 
C PLANE WAVE EXCITATION VECTOR ELEMENTS FOR 
WIRE AND 
C ATTACHMENT SEGMENT. INCIDENCE DIRECTION IS 
(THR,PHR). 
C - ONLY ONE INCIDENCE ANGLE PER CALL 
C -WIRE DOES NOT HAVE TO LIE IN THE Z=0 PLANE 
Cc 
COMPLEX 
U0, AA,BB,C,R(400),CEXP,EXP,FI1,Fl2,SI,D1,CMPLX 
DIMENSION 
ZH(400),NS(4),NW1(4),NW2(4),XH(400) YH(400) 
MP2=NP2-1 
MT2=NP2-2 
MT22=2°*MT2 
DO 5 L=1,NWIRES 
5 NS(LJ=NWe2(L)-NW1(L)+1 
Uo=(0.,0.) 
CC=COS(THR) 
SS=SIN(THR) 
CP=COS(PHR) 
SP=SIN(PHR) 
UPzSS°CP 
VP=Sorcr 
WP=CC 
DO 12 IP=1,MP2 
ll=IP+NP1 
l=|1+1 
ZS=.5*(ZH(1)+ZH(Il)) 
XS=.5*(XH(I)+XH(II)) 
YS=.5°(YH(I)+YH(II)) 
DX=XH(I)-XH(II) 
DY=YH(l)-YH(Il) 


D1=SQRT(DX**2+DY°**2) 
SU=DY/D1 
CU=DX/D1 
C FOR WIRES IN PERPENDICULAR TO Z SIN(V)=1 AND 
COS(V)=0 
SV=1.0 
Cv=0.0 
C 
C WIRE SEGMENT CALCULATIONS 
C 
A=UP*CU+VP*SU 
B=UP*XS+VP°YS+WP*ZS 
C=CMPLX(0.,A) 
EXP=CEXP(CMPLX(0.,B)) 
AA=CC*(CU*CP+SU*SP) 
BB=SU°CP-SP°CU 
PSl=D1°A/2. 
IF(PSILNE.0.) GO TO 60 
SINC=1. 
GO TO 61 
60 SINC=SIN(PSI)/PSI 
61 COSP=COS(PS)) 
FI1=SINC*D1°EXP/2. 
F1l2=(0.,0.) 
IF(ABS(A).LT.1.E-4) GO TO 62 
CSP=COSP-SINC 
IF(ABS(CSP).LT.1.E-4) GO TO 62 
FI2=EXP/C*CSP 
62 CONTINUE 
Sl=FI14+F 12 
Di=F11-Fl2 
Cc 
C R-WIRE-THETA 
Cc 
IF(IP.EQ.MP2) GO TO 10 
R(IP)=AA‘SI 
R(IP+MT2)=BB°SI 
10 CONTINUE 
C 
C R-WIRE-PHI 
C 
14 IF(IP.EQ.1) GO TO 12 
R(IP-1)=R(IP-1)+AA*DI 
R(IP-1+MT2)=R(IP-1+MT2)+BB"D! 
12 CONTINUE 
C 
C WRITE OVER DUMMY WIRE SEGMENTS 
o} 
IF(NWIRES.EQ.1) GO TO 210 
NSPTS=0 
DO 209 L=1,NWIRES-1 
NSPTS=NSPTS+NS(L) 
IN=NSPTS-2 
R(IN+1)=U0 
R(IN+2)=U0 
IN=IN+MT2 
R(IN+1)=U0 
R(IN+2)=U0 
c 209 CONTINUE 
210 RETURN 
END 
SUBROUTINE 
PLANES(M1,NPO,NP1,NP2,NT,RH.ZH,XT,AT,THR,R) 
Cc 


ooaonananngn0ono0on0agd0 a 
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C PLANE WAVE EXCITATION VECTOR FOR BOR 
ELEMENTS. 
C NO CHANGE FROM HARRINGTON'S ORIGINAL 
PROGRAM OTHER THAN 
Cc 
C °° ONLY ONE MODE PER CALL ** 
Cc 
COMPLEX 
R(1000),U,U1,UA,UB,FA(50),FB(50),F2A,F2B,F1A,F1B 
COMPLEX U2,U3,U4,U5,CMPLX 
DIMENSION RH(400),ZH(400) .XT(4),AT(4) 
DIMENSION R2(4),Z2(4),BJ(5000) 
W=(Ge 1) 
U1=3.141593°U°°M1 
NF=1 
M2=M1 
NP=NP0+NP14+NP2 
NRF=NPO+NP1 
MP=NP-1 
MT=MP- 1 
N=MT+MP 
N2=2*N 
CC=COS(THR) 
SS=SIN(THR) 
M3=M1+1 
M4=M2+3 
IF(M1.EQ.0) M3=2 
M5=M1+2 
M6=M2+2 
DO 12 IP=1,MP 
K2=IP 
l=IP+1 
DR=.5°(RH(I)-RH(IP)) 
DZ=.5°(ZH(I)-ZH(IP)}) 
D1=SQRT(DR*DR+D2Z°DZ) 
R1=.25°(RH(l)+RH(IP)) 
IF(ABS(R1).LT.1.E€-5) R1=1. 
Z21=.5°(ZH(I)+ZH(IP)) 
DR=.5°DR 
D2=DR/R1 
DO 13 L=1,NT 
R2(L)=R1+DR*XT(L) 
Z2(L)=Z1+DZ°XT(L) 
13 CONTINUE 
D3=DR*CC 
D4=-D2Z°SS 
DS=D1°CC 
DO 23 M=M3,M4 
FA(M)=0. 
FB(M)=0. 
23 CONTINUE 
DO 15 L=1,.NT 
M=Ss A2iL) 
IF(X.GT..5E-7) GO TO 19 
DO 20 M=M3,M4 
BJ(M)=0. 
20 CONTINUE 
Bu(Zi=T. 
Sei 
GO TO 18 
19 M=2.8°X+14.-2./X 
IF(X.LT..5) M=11.8+ALOG 10(X) 
IF(M.LT.M4) M=M4 
BJ(M)=0. 
JM=M-1 


BJU(JM)=1. 
DO 16 J=4,M 
J2=JM 
JM=JM-1 
J1=JM-1 
BU(JM)=J1/X°*BUJ(J2)-BU(JM+2) 
16 CONTINUE 
S=0. 
IF(M.LE.4) GO TO 24 
DO 17 J=4,M,2 
S=S+BJ(J) 
17 CONTINUE 
24 S=BJ(2)+2.°S 
18 ARG=Z2(L)°CC 
UA=AT(L)/S*CMPLX(COS(ARG),SIN(ARG)) 
UB=XT(L)*UA 
DO 25 M=M3,M4 
FA(M)=BJ(M)°UA+FA(M) 
FB(M)=BJ(M)*UB+FB(M) 
25 CONTINUE 
15 CONTINUE 
IF(M1.NE.0) GO TO 26 
FA(1)=- FA(3) 
FB(1)=- FB(3) 
26 UA=U1 
DO 27 M=M5,M6 
M7=N- 1 
M8=M+1 
F2A=UA*(FA(M8)+FA(M7)) 
FOB=UA*(FB(M8)+FB(M7)) 
UB=U"UA 
F1A=UB*(FA(M8)-FA(M7)) 
F1B=UB*(FB(M8)-FB(M7)) 
U4=D4°UA 
U2=D3°F1A+U4*FA(M) 
U3=D3°F1B+U4*FB(M) 
U4=DR*F2A 
U5=DR*F2B 
K1=K2-1 
K4=K1+N 
K5=K2+N 
R(K2+MT)=-D5*(F2A+D2*F2B) 
R(K5+MT)=D1°(F 1A+D2°F 1B) 
IF(IP.EQ.1) GO TO 21 
R(K1)=R(K1)+U2-U3 
R(K4)=R(K4)+U4-U5 
IF(IP.EQ.MP) GO TO 22 
21 R(K2)=U2+U3 
R(K5)=U4+U5 
22 K2=K2+N2 
UA=UB 
27 CONTINUE 
12 CONTINUE 
R(NPO- 1)=(0.,0.) 
R(NPO)=(0.,0.) 
R(NPO+MT)=(0.,0.) 
R(NRF-1)=(0.,0.) 
R(NRF)=(0.,0.) 
R(NRF+MT)=(0.,0.) 
99 RETURN 
END 
SUBROUTINE 
ZMATSS(M1,M2,NP0,NP1,NP2,NPHI,NT,RH.ZH,X,A,XT,A 
TZ 
Cc 
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C ****ee** THREE DETACHED SURFACES °*********** 
C SURFACE-SURFACE IMPEDANCE ELEMENTS. 
REMAINS UNCHANGED FROM HARRINGTON 
C EXCEPT MULTIPLE SURFACES ARE PERMITTED. 
THE FIRST SURFACE HAS NPO 
C POINTS, THE SECOND NP1 POINTS AND THE THIRD 
NP2. 
C NPO PARABOLOID PTS; NP1 HYPERBOLOID PTS; 
NP2 CAVITY PTS. 
Cc 
COMPLEX 
2(50000),U1,U2,U3,U4,U5,U6,U7,U8,U9,UA,UB,G4A(4),G5 
A(4) 
COMPLEX 
G6A(4),G4B(4),G5B(4),G6B(4),H4A,H5A,H6A,H4B,H5B 
COMPLEX CMPLX,H6B,UC,UD,GA(400),GB(400) 
DIMENSION 
RH(400),ZH(400),X(400),A(400),AT(4),RS(400),ZS(400) 
DIMENSION 
D(400),DR(400),D2Z(400),DM(400) ,C2(400),C3(400),R2(4) 
DIMENSION 
C4(400),C5(400),C6(400),27(4),R7(4),Z8(4),R8(4) 
DIMENSION XT(4),22(4) 
Cil=2: 
Cr= 
NP=NP0+NP14+NP2 
NRF=NPO+NP1 
DO 10 l=2,NP 
I2=1-1 
RS(I2)=.5*(RH(I)}+RH(12)) 
ZS(I2)=.5*(ZH(1)+ZH(12)) 
D1=.5°(RH(I)-RH(I2)) 
D2=.5°(ZH(I)-ZH(I2)) 
D(I2)=SQRT(D1°D1+D2°D2) 
DR(I2)=D1 
D2Z(I2)=D2 
DM(I2)=D(I2)/RS(12) 
10 CONTINUE 
M3=M2-M1+1 
M4=M1-1 
PIl2=1.570796 
DO 11 K=1,NPHI 
PH=PI2°(X(K)+1.) 
C2(K)=PH*PH 
SN=SIN(.5°PH) 
C3(K)=4.°SN°SN 
A1=PI2*A(K) 
D4=.5°A1°C3(K) 
D5=A1*COS(PH) 
D6=A1*SIN(PH) 
M5=K 
DO 29 M=1,M3 
PHM=(M4+M)*°PH 
A2=COS(PHM) 
C4(M5)=D4"A2 
C5(M5)=D5°A2 
C6(M5)=D6*SIN(PHM) 
M5=M5+NPHI 
29 CONTINUE 
11 CONTINUE 
MP=NP-1 
MT=MP-1 
N=MT+MP 
N2N=MT*N 
N2=N*N 
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Uq=(0.,.5) 

U2=(0.,2.) 

JN=-1-N 

DO 15 JOQ=1,MP 
KQ=2 

IF(JO.EQ.1) KQ=1 
IF(JQ.EQ.MP) KQ=3 
R1=RS(JQ) 
Z21=ZS(JQ) 
D1=D(JQ) 
D2=DR(JQ) 
D3=DZ(JQ) 
D4=D2/R1 
D5=DM(JQ) 
SV=D2/D1 
CV=D3/01 
T6=CT‘D1 
T62=T6+D1 
T622T62°T62 
R6=CP*R1 
R62=R6°R6 

DO 12 L=1,NT 
R2(L)=R1+D2°XT(L) 
Z2(L)=Z1+D3°XT(L) 
12 CONTINUE 
U3=D2°U1 
U4=D3°U1 

DO 16 IP=1,MP 
R3=RS(IP) 
Zs=2 3 IP) 
R4=R1-R3 
Z24=Z21-Z3 
FM=R4°SV+Z4°CV 
PHM=ABS(FM) 
PH=ABS(R4°CV-2Z4°SV) 
D6=PH 
IF(PHM.LE.D1) GO TO 26 
D6=PHM-D1 
D6=SORT(D6*D6+PH*PH) 
26 IF(IP.EQ.JQ) GO TO 27 
Ke 

IF(T6.GT.D6) KP=2 
IF(R6.GT.D6) KP=3 
GO TO 28 


27 KP=4 

28 GO TO (41,42,41,42),KP 

42 DO 40 L=1,NT 
D7=R2(L)-R3 
D8=Z2(L)-Z3 
Z7(L)=D7°D7+D8*D8 
R7(L)=R3*R2(L) 
Z8(L)=.25°Z7(L) 
R8(L)=.25°R7(L) 

40 CONTINUE 
Z4=R4°R44+2Z4°24 
R4=R3°R1 
R5=.5°R3°SV 
DO 33 K=1,NPHI 
A1=C3(K) 

RR=Z24+R4°A1 

UA=0. 

UB=0. 

IF(RR.LT.T62) GO TO 34 
DO 35 L=1,NT 
R=SQRT(Z7(L)+R7(L)°A1) 
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SN=-SIN(R) 
CS=COS(R) 
UC=AT(L)/R°CMPLX(CS,SN) 
UA=UA+UC 
UB=XT(L)*UC+UB 
35 CONTINUE 
GOTO 36 
34 DO 37 L=1,NT 
R=SQRT(Z8(L)+R8(L)*A1) 
SN=-SIN(R) 
CS=COS(R) 
UC=AT(L)/R*SN*CMPLX(-SN,CS) 
UA=UA+UC 
UB=XT(L)*UC+UB 
37 CONTINUE 
A2=FM+R5°A1 
D9=RR-A2*A2 
R=ABS(A2) 
D7=R-D1 
D8=R+D1 
D6=SORT(D8*D8+D9) 
R=SQRT(D7*D7+D9) 
IF(D7.GE.0.) GO TO 38 
A1=ALOG((D8+D6)*(-D7+R)/D9)/D1 
GO TO 39 
38 A1=ALOG((D8+D6)/(D7+R))/D1 
39 UA=A14+UA 
UB=A2°(4./(D6+R)-A1)/D14+UB 
36 GA(K)=UA 
GB(K)=UB 
33 CONTINUE 
K1=0 
DO 45 M=1,M3 
H4A=0. 
HS5A=0. 
H6A=0. 
H4B=0. 
H5B=0. 
H6B=0. 
DO 46 K=1,NPHI 
K1=K1+1 
D6=C4(K1) 
D7=C5(K1) 
D8=C6(K1) 


UA=GA(K) 

UB=GB(K) 

H4A=D6*UA+H4A 
H5A=D7*UA+H5A 
H6A=D8*UA+H6A 
H4B=D6*UB+H4B 
H5B=D7*UB+H5B 
H6B=D8*UB+H6B 


46 CONTINUE 


G4A(M)=H4A 
G5A(M)=H5A 
G6A(M)=HEA 
G4B(M)=H4B 
G5B(M)=H5B 
G6B(M)=H6B 


45 CONTINUE 


IF(KP.NE.4) GO TO 47 
=D1/(P12°R1) 

D6=2./01 

Ds=0. 

DO 63 K=1,NPHI 


A1=R4°C2(K) 
R=R4°C3(K) 
IF(R.LT.T62) GO TO 64 
D7=0; 
DO 65 L=1,NT 
D7=D7+AT(L)/SQRT(Z7(L)+A1) 
65 CONTINUE 
GO TO 66 
64 A1=A2/(X(K)}+1.) 
D7=D6*ALOG(A14+SQRT(1.4+A1°A1)) 
66 D8=D8+A(K)°D7 
63 CONTINUE 
A1l=.5°A2 
A2=1./A1 
D8=-P12°D8+2./R1*(BLOG(A2)+A2°BLOG(A1)) 
DO 67 M=1,M3 
G5A(M)=D8+G5A(M) 
67 CONTINUE 
GO TO 47 
41 DO 25 M=1,M3 
G4A(M)=0. 
G5A(M)=0. 
G6A(M)=0. 
G4B(M)=0. 
G5B(M)=0. 
G6B(M)=0. 
25 CONTINUE 
DO 13 L=1,NT 
A1=R2(L) 
R4=A1-R3 
24=Z2(L)-23 
Z24=R4°R4+24°24 
R4=R3°A1 
DO 17 K=1,NPHI 
R=SORT(Z4+R4°C3(K)) 
SN=-SIN(R) 
CS=COS(R) 
GA(K)=CMPLX(CS,SN)/R 
17 CONTINUE 
DS=0, 
IF(R62.LE.24) GO TO 51 
DO 62 K=1,NPHI 
D6=D6+A(K)/SQRT(Z4+R4°C2(K)) 
62 CONTINUE 
Z24=3.141593/SORT(Z4/R4) 
D6=-P12°D6+ALOG(Z4+SQRT(1.+24°Z4))/SQRT(R4) 
51 A1=AT(L) 
A2=XT(L)*A1 
K1=0 
DO 30 M=1,M3 
W5=0; 
U6=0. 
U7=0. 
DO 32 K=1,NPHI 
UA=GA(K) 
K1=K1+1 
U5=C4(K1)°UA+U5 
U6=C5(K1)°UA+U6 
U7=C6(K1)*UA+U7 
32 CONTINUE 
U6=D6+U6 
G4A(M)=A1°U5+G4A(M) 
G5A(M)=A1°U6+G5A(M) 
G6A(M)=A1°U7+G6A(M) 
G4B(M)=A2°U5+G4B(M) 
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G5B(M)=A2*U6+G5B(M) 
G6B(M)=A2°U7+G6B(M) 
30 CONTINUE 
13 CONTINUE 
47 A1=DR(IP) 
UA=A1°U3 
UB=D2(IP)*U4 
A2=D(IP) 
D6=-A2*D2 
D7=D1°A1 
D8=D1°A2 
JM=JN 
DO 31 M=1,M3 
FM=M4+M 
A1=FM*DM(IP) 
H5A=G5A(M) 
H5B=G5B(M) 
H4A=G4A(M)+H5A 
H4B=G4B(M)+H5B 
H6A=G6A(M) 
H6B=G6B(M) 
U7=UA°H5A+UB*H4A 
U8=UA°H5B+UB°H4B 
U5=U7-U8 
U6=U7+U8 
U7=-U1°H4A 
U8=D6*H6A 
U9=D6°H6B-A1°H4A 
UC=D7*(H6A+D4°H6B) 
UD=FM*D5*H4A 
K1=IP+JM 
K2=K1+1 
K3=K1+N 
K4=K2+N 
K5=K2+MT 
K6=K44+MT 
K7=K3+N2N 
K8=K44N2N 
GO TO (18,20,19),KQ 
18 Z(K6)=U8+U9 
IF(IP.EQ.1) GO TO 21 
Z(K3)=Z(K3)+U6-U7 
2(K7)=Z(K7)+UC-UD 
IF(IP.EQ.MP) GO TO 22 
21 Z(K4)=U6+U7 
Z(K8)=UC+UD 
GO TO 22 
19 Z(K5)=Z2(K5)+U8-U9 
IF(IP.EQ.1) GO TO 23 
Z2(K1)=Z(K1)+U5+U7 
2(K7)=Z(K7)+UC-UD 
IF(IP.EQ.MP) GO TO 22 
23 Z(K2)=Z(K2)+U5-U7 
Z(K8)=UC+UD 
GO TO 22 
20 2(K5)=Z(K5)+U8-U9 
Z(K6)=U8+U9 
IF(IP.EQ.1) GO TO 24 
Z(K1)=Z(K1)4+U5+U7 
Z(K3)=Z(K3)+U6-U7 
Z(K7)=Z(K7)+UC-UD 
IF(IP.EQ.MP) GO TO 22 
24 Z2(K2)=Z(K2)+U5-U7 
Z2(K4)=U6+U7 
2(K8)=UC+UD 


22 Z(K8+MT )=U2*(D8*(H5A+D4°H5B)-A1°UD) 
JM=JM4+N2 

31 CONTINUE 

16 CONTINUE 
JN=JN+N 
15 CONTINUE 
Cc 
C THREE MULTIPLE BODIES USING THE SIMPLIFIED 
APPROACH 
C NULL OUT ROWS AND COLS FOR THE FIRST 
DUMMY SEGMENT 
C 
DO 100 LSS=1,N 
ES=EsS-1 
Z(LS*N+NP0O)=(0.,0.) 
Z(LS*N+NPO-1)=(0.,0.) 

100 Z(LS*N4+NP0+MT)=(0.,0.) 

BS 101 ES=1,N 
Z((NPO-2)°N+LS)=(0.,0.) 
Z((NPO-1)°N+LS)=(0.,0.) 

101 Z((NPO-1+MT)*N+LS)=(0.,0.) 
Z((NPO-2)*N+NPO- 1)=(1.,0.) 
Z((NPO-1)*N+NP0)=(1.,0.) 
Z((MT+NPO-1)"N+NPO+MT)=(1.,0.) 

Se 
C NULL OUT ROWS AND COLS FOR THE SECOND 
DUMMY SEGMENT 
Cc 
DO 200 LSS=1,N 
LS=LSS-1 
Z(LS*N+NRF)=(0.,0.) 
Z(LS*N+NRF- 1)=(0.,0.) 

200 Z(LS*N4+NRF+MT)=(0.,0.) 

BO 201 LS=1,N 
Z((NRF-2)*N+LS)=(0.,0.) 
Z((NRF-1)°N+LS)=(0.,0.) 

201 Z((NRF-14+MT)*N+LS)=(0.,0.) 
Z((NRF-2)°N+NRF-1)=(1.,0.) 
Z((NRF-1)°N+NRF)=(1.,0.) 
Z((MT+NRF- 1)°N+#+NRF+MT)=(1.,0.) 

999 RETURN 
END 
SUBROUTINE ZASMBO(NP1,NP2,MODES,Z,ZWW) 

C 

C ASSEMBLES MODE INDEPENDENT BLOCKS OF Z 

Cc 
COMPLEX 2Z(500000) ZWW( 10000) 
MT1=NP1-2 
MP1=NP1-1 
N=MT14+MP1 
MT2=NP2-2 
NBLOCK=2*MODES +1 
NS=NBLOCK‘'N 
NR=NS4+MT2 
NS2=NS‘NR 
NWWE=NS2+NS 
DO 10 iC=1,MT2 
iC 1=IC-1 
DO 10 IR=1,MT2 
K=NWW+IC 1°NR+IR 
IWeIC1°MT2+IR 

10 2Z(K)=ZWW(IW) 
RETURN 
END 
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SUBROUTINE 
ZASMBN(NP1,NP2,MODES,.NM,Z.ZSS.ZSW) 
C 
C ASSEMBLES MODE DEPENDENT BLOCKS IN Z. 
MODES=TL # OF MODES; 

C NM IS THE MODE # OF THE BLOCKS ZSS.ZSW.,ZSJ. 
(IF NM<0 THEN 

C ZSS IS OMITTED - ONLY FILLED FOR POSITIVE NM) 
Cc 

COMPLEX ZSS(50000),2(500000), ZSW(50000) 

MT1=NP 1-2 

MP 1=NP1-1 

MT2=NP2-2 

N=MP 1+MT 1 

NBLOCK=2*MODES+1 

NS=NBLOCK’*N 

NR=NS+MT2 

MT=MT14+MT2 

MTN=MT°N 

IF(NM.LT.0) GO TO 225 
C 
C ZSS BLOCKS FOR MODE NM (+ AND - FILLED IN 1 
LOOP) 

Cc 

NTOP1=MODES-NM 

NTOP2=NBLOCK-(NTOP1+1) 

NSS2=NTOP 1*(N4+NR°‘N) 

NSS1=NTOP2*(N+NR‘N) 

DO 200 [C=1,MP1 

IC 1=IC-1 

DO 200 IR=1,MP1 

IF(IC.GT.MT1) GO TO 100 

IF(IR.GT.MT1) GO TO 50 

IO=IC1°N+IR 

ITT=IC 1°NR+IR 

Z(ITT+NSS1)=ZSS(I0) 

IF(NM.NE.O) Z(ITT+NSS2)=ZSS(I0) 
50 l0=MT1+IC1°N+IR 

IPT=IC1°NR4+MT141R 

Z(IPT+NSS1)= ZSS(I0) 

IF(NM.NE.O) Z(IPT4+#NSS2)=-ZSS(10) 
100 IF(IR.GT.MT1) GO TO 150 

IO=N*MT1+IC1°N+IR 

ITP=MT 1°NR+IC1°NR+IR 

Z(ITP+NSS1)= ZSS(I0) 

IF(NM.NE.0) Z(ITP+NSS2)=-ZSS(I0) 
150 10=N*MT14+MT1+IC1°N+IR 

IPP=IC1*°NR+MT14MT1°NR+IR 

Z(IPP+NSS1)=ZSS(I0) 

IF(NM.NE.O) Z(IPP+NSS2)=ZSS(I0) 
200 CONTINUE 
225 CONTINUE 
C 
C ZSW AND ZWS BLOCKS 
cS 

NTOP1=MODES+NM 

NTOP2=NBLOCK-(NTOP1+1) 

NSW=NTOP1°N4+NS*NR 

NWS=NTOP2°NR'N4+NS 

DO 250 IC=1,MT2 

IC1=IC-1 

DO 250 IR=1,N 

ISW=IR+IC1°NR+NSW 

IWS=(IR-1)*NR+IC+NWS 

lO=IC1°N+IR 


Z(ISW)=ZSW(I0) 
Z(IWS)=ZSW(I0) 
250 CONTINUE 
RETURN 
END 
SUBROUTINE 
ZMATSW(NWIRES,NW1,NW2,NP11,NP12,NP 13, 
* XH,YH,RH,ZH,NT,XT,AT,NG,XG,AG,MODE,A,Z) 
Cc 
C >>>>>>>>>>> revised version <<<<<<<<<< 
Cc 
C CALCULATION OF THE SURFACE-WIRE MATRIX 
ELEMENTS. ATTACHMENT 
C POINT NOT AT A SURFACE EDGE, BUT CAN BE ON Z 
AXIS (IC#1). SURFACE 
C START INDEX IS 1. TWO BORS OF LENGTHS NP11, 
NP12 AND NP 13. WIRE 
C START AND STOP POINTS ARE NW1 AND NW2. 
Cc 
COMPLEX 
CEXP,EXP,2Z(50000),G5,G6,G7,CONH,F5,F6,F7 
COMPLEX 
U,U0,PSI,T0,T1,12,73,ST1,51T2,ST3,CON,.SP1,SP2.SP3 
COMPLEX U5,U6,U7,U8,U9,CMPLX,T4,T5,T6,T7 
DIMENSION 
RS1(400),RH(400),ZH(400),XG(400), AG(400),XT(4),AT(4) 
DIMENSION 
ZS 1(400),01(400),0(400),S(400),C(400),RS(400) 
DIMENSION 
ZS(400),S1(400),C 1(400), NW1(4), NW2(4),NS(4) 
DIMENSION 
XH(400),YH(400),XS1(400),¥S1(400),CU(400),SU(400) 
Pl=3.14159 
FiZz=2.r) 
NP1=NP11+NP12+NP13 
NRF=NP11+NP12 
DO 1L=1,NWIRES 
1 NS(L)=NW2(L)-NW1(L)+1 
Ud=(0.,0.) 
WStOei) 
Gg 
C DEFINE SURFACE GEOMETRY TERMS 
Cc 
MP=NP 1-1 
MT=MP-1 
NROW=MP+MT 
DO 5 l=2,NP 1 
I2=I-1 
RS(1l2)=.5°(RH(I)+RH(12)) 
ZS9(12)=.5*°(ZH(l)+ZH(I2)) 
DR=RH(1)-RH(I2) 
DZ=ZH(1)-ZH(I2) 
D0(I2)=SQRT(DR**2+DZ**2) 
$(I2)=DR/D(I2) 
5 C(l2)=DZ/D(I2) 
C 
C DEFINE WIRE GEOMETRY TERMS 
Cc 
NW=NW2(NWIRES)-NW1(1) 
NWP=NW+1 
MT2=NWP-2 
DO 10 N=2,NWP 
NO=N-1 
I=NP14+N 
I2=l-1 
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RS1(NQ)=.5*(RH(I)+RH(I2)) 
Z2S1(NO)=.5*(ZH(1)+ZH(12)) 
XS 1(NOQ)=.5°(XH(l)+XH(I2)) 
YS1(NO)=.5*(YH(1)+YH(I2)) 
DX=XH(I)-XH(I2) 
DY=YH(I)-YH(I2) 
D1(NQ)=SQRT(DX**2+DY**2) 
C FOR WIRES PARALLEL TO THE XY PLANE SIN(V)=1 
AND COS(V)=0 
$1(NQ)=1. 
C1(NQ)=0. 
CU(NO0)=DX/D1(NO) 
SU(N0)=DY/D1(NO) 
10 CONTINUE 
CON=SUPI2/2. 
C INTEGRATING FROM 0 TO 2°PI 
O1=PIl 
O2=P1 
JN=-1-NROW 
C 
C SEGMENT LOOPS: g=WIRE, p=SURFACE 
C 
DO 100 JQ=1,NW 
KQ=2 
IF{JG EO 1) KQ=1 
IF(JQ.EQ.NW) KQ=3 
H1=01(JQ)/2. 
ZW=ZS1(JQ) 
DO 22 IP=1,MP 
Si l=uo 
ST2=U0 
ST3=U0 
SP1=U0 
SP2=U0 
Sea=u0 
AA=S(IP)*CU(JQ) 
BB=S(IP)*SU(JQ) 
RD=RS(IP)*D1(JQ) 
DD=D(IP)*D1(JQ) 
CONH=CON’D(IP) 
Cc 
C FIRST TERM (H-INTEGRATION) 
Cc 
DO 30 l=1,NT 
H=H 1*XT (I) 
HHD=2.°H/D1(JQ) 
XWeXS1(JQ)+H*CU(JQ) 
YWe=YS1(JQ)+H*SU(JQ) 
C 
C PHI INTEGRATIONS FOR G-HAT FUNCTIONS 
Cc 
G7=U0 
G5=U0 
G6=U0 
DO 32 K=1,NG 
PH=01*XG(K)+02 
PHM=PH*FLOAT(MODE) 
XS=RS(IP)*COS(PH) 
YS=RS(IP)*SIN(PH) 
EXP=CEXP(CMPLX(0.,-PHM)) 
RP=SQRT((XS-XW)**2+(YS-YW)**24(ZS(IP)-ZW) **2) 
PSI=CEXP(CMPLX(0.,-RP))/RP 
F5=EXP*COS(PH)*PSI| 
F6=EXP*SIN(PH)‘PSI 
F7=EXP*PSI 


G7=G7+F7*AG(K) 
G5=G5+F5*AG(K) 
32 G6=G6+F6*AG(K) 
G7=01°G7 
G5=01°G5 
G6=01°G6 
T0=G5°AA/4. 
T1=TO*HHD 
T2=G6°BB/4. 
Té=T2°HHD 
T4=G5°SU(JQ)/2. 
T5=T4*HHD 
T6=-G6*CU(JQ)/2. 
T7=T6*HHD 
ST 1=ST1+AT(I)*(T0+T2) 
ST2=ST2+AT(I)*(T1+T3) 
ST3=ST3+AT(I)}°G7/DD 
SP1=SP1+AT(I)°(T4+T6) 
SP2=SP2+AT(I)*(15+T7) 
SP3=SP3+AT(I)°U°G7°FLOAT(MODE)/RD 
30 CONTINUE 
Cc 
C ZSW-t TERMS (K1,K2,K3,K4) 
Cc 
U5=(ST1-ST2)*CONH°H1 
U6=(ST1+ST2)°CONH*H1 
U7=-ST3°CONH*H1 
Cc 
C ZSW-p TERMS (K5,K6) 
Cc 
U8=SP 1°CONH‘H1 
U9=(SP2+SP3)*CONH’H1 
K1i=IP4+JN 
K2=K1+1 
K3=K1+NROW 
K4=K2+NROW 
K5=K2+MT 
K6=K4+MT 
GO TO (18,20,19),KQ 
18 Z(K6)=U8+U9 
IF(IP.EQ.1) GO TO 21 
Z(K3)=Z(K3)+U6-U7 
IF(IP.EQ.MP) GO TO 22 
21 Z(K4)=U6+U7 
GO TO 22 
19 Z(K5)=Z(K5)+U8-U9 
IF(IP.EO.1) GO TO 23 
Z(K1)=Z2(K1)+U5+U7 
IF(IP.EQ.MP) GO TO 22 
23 Z(K2)=Z2(K2)+U5-U7 
GO TO 22 
20 Z(K5)=Z(K5)+U8-U9 
Z(K6)=U8+U9 
IF(IP.EQ.1) GO TO 24 
Z(K1)=Z(K1)+U5+U7 
Z(K3)=2(K3)+U6-U7 
IF(IP.EQ.MP) GO TO 22 
24 Z(K2)=Z2(K2)+U5-U7 
Z(K4)=U64+U7 
22 CONTINUE 
JN=JN+NROW 
100 CONTINUE 
Cc 
C FOR MULTIPLE BORS WRITE OVER DUMMY 
SEGMENT ROWS 
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C WRITE OVER ROWS AND COLS FOR FIRST DUMMY 
SEG. 
Cc 
DO 101 LSS=1,MT2 
LS=LSS-1 
2(LS°NROW+NP11)=(0.,0.) 
Z(LS*NROW+NP11-1)=(0.,0.) 
101 Z(LS°NROW+NP114MT)=(0.,0.) 
é 
C WRITE OVER ROWS AND COLS FOR SECOND 
DUMMY SEG. 
C 
DO 201 LSS=1,MT2 
LS=LSS-1 
Z(LS*°NROW+NRF)=(0.,0.) 
Z(LS*°NROW+NRF- 1)=(0.,0.) 
201 Z(LS°NROW+NRF+MT)=(0.,0.) 
Cc 
C WRITE OVER DUMMY SEGMENTS IN THE CASE OF 
MULTIPE WIRES 
Cc 
IF(NWIRES.EQ.1) GO TO 292 
NSPTS=0 
DO 291 L=1,NWIRES-1 
NSPTS=NSPTS+NS(L) 
DO 291 l=1,NROW 
c 
C SET COLS TO ZERO FOR STRIPS 
c 
IN=(NSPTS-2)°NROW 
Z(IN+1)=U0 
Z(IN+1+NROW)=U0 
291 CONTINUE 
292 RETURN 
END 
SUBROUTINE 
ZMATWW(NWIRES,NW1,NW2,XH,YH,RH.ZH,NT,XT,AT,A, 
2) 
C 
C *** MODS FOR DIPOLE -- RECTANGULAR 
COORDINATES ARE SENT OVER *** 
C COMPUTE MATRIX ELEMENTS FOR MULTIPLE 
WIRES USING GLISSON METHOD 
C NWIRES= # OF WIRES; NW1(I), NW2(1) ARE START 
AND STOP POINTS. 
C &&&&& WIRES NEED NOT BE IN THE Z=0 PLANE 
&8888 
c 
COMPLEX CEXP,Z(10000),CON,CMPLX 
COMPLEX U,U0,PSI,SUM1,SUM2,SUM3,U5,U6,U7 
DIMENSION 
RH(400),ZH(400),XT(4),AT(4),XH(400), YH(400) 
DIMENSION D1(400),$1(400),C1(400) 
DIMENSION 
NW1(4),NW2(4),NS(4),XS1(400), YS 1(400),ZS1(400) 
DIMENSION CU(400),SU(400) 
Pl=3.141592 
Pi2=2.*P| 
PD2=P1/2. 
U0=(0.,0.) 
U=(0.,1.) 
c 
C DEFINE GEOMETRY TERMS FOR THE WIRE. 
XH,YH,RH,ZH ARE ALL KNOWNS. 
CG 


DO 5 L=1,NWIRES 

5 NS(LJ=fENW2(L)-NW1(L)+1 
NS1=NW2(NWIRES)-NW1(1) 
NPS=NS1+1 
DO 10 N=2,NPS 
NO=N-1 
l=NW1(1)+N-1 
I2=!-1 
XS1(NO)=.5°*(XH(I)+XH(12)) 
YS1(NO)=.5°(YH(I)+YH(I2)) 
ZS1(NO)=.5°(ZH(1}+ZH(I2)) 
DX=XH(I)-XH(I2) 
DY=YH(I)- YH(12) 
D1(N0)=SORT(DX°**2+DY*°2) 
CU(N0O)=DX/D1(NO) 
SU(NO)=DY/D 1(NO) 
$1(NO)=1.0 
C1(N0O)=0.0 

10 CONTINUE 
NROWENS1-1 
JN=-1-NROW 
DO 500 JO=1,NS1 
KO=2 
IF(JO.EO.1) KO=1 
IF(JO.EO.NS1) KO=3 


01=D1(JO)/2. 
DO 22 IP=1,NS1 
LO=0 


IF(IP.EO.JO) LO=1 
CON=D1(IP)/4./PI°U 
SUM1=U0 
SUM2=U0 
SUM3=U0 
AA=CU(IP)*CU(JO)+SU(IP)*SU(JO) 
C 
CH INTEGRATION -- SUBTRACT OUT SINGULARITY IF 
JO=IP 
C AND SAME SEGMENT (N11=N21): DESIGNATED LO=1 
Cc 
DO 100 J=1,NT 
H=01°XT(l) 
IF(LO.NE.1) GO TO 40 
RP=SORT(H**2+A°A) 
GO TO 50 
40 XP=XS1(JO)+H"*CU(JO)*S 1(JO) 
YP=YS 1(JO)+H*SU(JO)*S 1(JO) 


RP=SORT((XS1(IP)-XP)°°2+(YS1(IP)-YP)°°2+(ZS1(JO)-ZS 

1(IP))°*2) 

50 PSI=CEXP(CMPLX(0.,-RP)) 

c IF(LO.NE.1) GO TO 60 

¢ PSIl=PSI-(1.,0.) 

60 PSI=PSI/RP 
SUM1=SUM1+PSI°AA*AT(I)/4. 
SUM2=SUM2+PSI*AA*AT(I)/2./D1(JO)*H 
SUM3=SUM3+PSI*°AT(I)/D1(IP}/D1(JO) 

100 CONTINUE 
X120. 

2-0: 

X3=0. 

IF(LO.NE.1) GO TO 200 
DO=D1(JO)/2. 
SDO=SORT(DO*"2+A**2) 
XI0=ALOG(SDO+D0O)-ALOG(SDO-DO) 
X1=X10°AA/4. 


oo ona 0 
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c¢ X3=X10/D1(IP)/D1(UO) 
¢ 200 CONTINUE 
SUM1=(SUM1*°O1+X1)°CON 
SUM2=(SUM2°O1+X2)"CON 
SUM3=(SUM3°01+X3)"CON 
U5=SUM1-SUM2 
U6=SUM1+SUM2 
U7=-SUM3 
K1sIP+JN 
K2=K1+1 
K3=K1+NROW 
K4=K2+NROW 
GO TO (18,20,19),KO 
18 CONTINUE 
IF(IP.EO.1) GO TO 21 
Z(K3)=Z(K3)+U6-U7 
IF(IP.EQ.NS1) GO TO 22 
21 Z(K4)=U6+U7 
GOTO 22 
19 CONTINUE 
IF(IP.EO.1) GO TO 23 
Z(K1)=Z(K 1)+U54+U7 
IF(IP.EO.NS1) GO TO 22 
23 Z(K2)=Z(K2)+U5-U7 
GO TO 22 
20 CONTINUE 
IF(IP.EO.1) GO TO 24 
Z(K1)=Z(K1)+U54+U7 
Z(K3)=Z(K3)+U6-U7 
IF(IP.EO.NS1) GO TO 22 
24 Z(K2)=Z(K2)+U5-U7 
Z(K4)=U6+U7 
22 CONTINUE 
JN=JIN+NROW 
500 CONTINUE 
IF(NWIRES.EO.1) GO TO 297 
NSPTS=0 
DO 291 L=1,NWIRES-1 
NSPTS=NSPTS+NS(L) 
DO 291 l=1.NROW 
Cc 
C USE THE SIMPLIFIED METHOD OF HANDLING 
MULTIPLE WIRES 
C SET COLS TO ZERO FOR STRIPS 


Cc 
IN=(NSPTS-2)*NROW 
Z(IN+I}=U0 
Z(IN+I4+NROW)=U0 

C 


C SET ROWS TO ZERO FOR STRIPS 

e 
IN-NSPTS-2+(I-1)"NROW 
Z(IN+1)=U0 
Z(IN+2)=U0 

291 CONTINUE 

G 

C SET DIAGONALS TO 1 FOR STRIPS 

Cc 
NSPTS=0 
DO 296 L=1,NWIRES- 1 
NSPTS=NSPTS+NS(L) 
l0=NSPTS-2 
11=I0°NROW+NSPTS-2 
Z(11+1)=(1.,0.) 
Z(114+NROW+2)=(1.,0.) 


296 CONTINUE 
297 RETURN 
END 
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